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We discuss non-Gaussian effects associated with the local large-scale structure contributions to the 
Cosmic Microwave Background (CMB) anisotropics through the thermal Sunyaev-Zel'dovich (SZ) 
effect. The non-Gaussianities associated with the SZ effect arise from the existence of a significant 
four-point correlation function in large scale pressure fluctuations. Using the pressure trispectrum 
calculated under the recently popular halo model, we discuss the full covariance of the SZ thermal 
power spectrum. We use this full covariance matrix to study the astrophysical uses of the SZ effect 
and discuss the extent to which gas properties can be derived from a measurement of the SZ power 
spectrum. With the SZ thermal effect separated in temperature fluctuations using its frequency in- 
formation, the kinetic SZ effect, also known as the Ostriker-Vishniac effect, is expected to dominate 
the CMB temperature fluctuations at small angular scales. This effect arises from the baryon modu- 
lation of the first order Doppler effect resulting from the relative motion of scatterers. The presence 
of the SZ kinetic effect can be determined through a cross-correlation between the SZ thermal and 
a CMB map at small scales. Since the SZ kinetic effect is second order, however, contributions to 
such a cross-correlation arise to the lower order in the form of a three-point correlation function, 
or a bispectrum in Fourier space. We suggest an additional statistic that can be used to study 
the correlation between pressure traced by the SZ thermal effect and the baryons traced by the SZ 
kinetic effect involving the cross-power spectrum constructed through squared temperatures instead 
of the usual temperature itself. Through a signal-to-noise calculation, we show that future small 
angular scale multi-frequency CMB experiments, sensitive to multipoles of a few thousand, will be 
able to measure the cross-correlation of SZ thermal and SZ kinetic effect through a temperature 
squared power spectrum. 
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I. INTRODUCTION 

In recent years, motivated by the advances in the ex- 
perimental front, increasing attention has been given to 
theoretical details related to secondary anisotropics in 
the cosmic microwave background (CMB) resulting from 
the low redshift large scale structure. In addition to the 
direct detection, some of these effects have also been 
studied in detail due to their non-linear behavior lead- 
ing to higher order correlations in the CMB temperature 
fluctuations |l|,||. In general, the secondary anisotropies 
resulting from low redshifts {z ^ few tens) can be divided 
into those that result from gravitational effects and those 
resulting from scattering via free electrons in the rcion- 
ized epoch. 

Through the imparted differential gravitational red- 
shift of CMB photons, when traversing time-evolving 
potential fluctuations, the integrated Sachs- Wolfe (ISW; 
1^) contributes at large angular scales. This effect is 
an important contributor at z ^ 1 in a ACDM uni- 
verse. Through angular deflections across the sky via 
foreground gravitational lensing potentials, weak grav- 
itational lensing modifies the CMB power spectrum at 
large angular scales while creating power primarily at 



small angular scales (e.g., [Q-^ and references therein). 
This lensing effect also correlates with potentials traced 
by other secondary effects and produces a bispectrum in 
CMB temperature data (e.g., ||l|,||). 

In terms of the scattering effects, the inverse- Compton 
scattering of CMB photons via hot electrons, the so- 
called Sunyaev-Zel'dovich (SZ; |^) effect, dominates the 
small angular scale signal. The SZ effect has now been 
directly imaged towards massive galaxy clusters (e.g., 
|]lO| , p| ), where the temperature of the scattering medium 
can reach as high as 10 keV, producing temperature 
changes in the CMB of order 1 mK at Rayleigh-Jeans 
wavelengths. Given that the SZ effect also bears a spec- 
tral signature that differs from true temperature fluctua- 
tions, SZ contribution can be separated in multifrequency 
data, thereby suggesting a means of studying its proper- 
ties. As discussed in detail in |l^, a multi-frequency 
approach can be applied to current Boomerang [T^ , and 
upcoming MAP0 and Planck surveyor^ missions. 

The improving capabilities of hydrodynamical simula- 
tions and analytical models involving the large scale gas 
distribution have now provided predictions for the SZ ef- 
fect, mainly the SZ power spectrum and the expected 
number counts of SZ halos (see, |T2|jr^ -p^ among vari- 



*http:/ /map. nasa.gsfc.gov 

■|^http://astro.estec.esa.nl/Planck/; also, ESA D/SCI(6)3. 



1 



ous other studies). A wide-field SZ image, allowing de- 
tailed statistical studies such as the angular power spec- 
trum, is yet to be produced, though several experimen- 
tal attempts are currently in progress for this purpose. 
These experiments include the interferometric survey at 
the combined BIMA/OVRO array (CARMA; John Carl- 
strom, private communication) , the MINT interferometer 
(Lyman Page, private communication), and the BOLO- 
CAM array on the Caltech Submillimeter Observatory 
(Andrew Lange, private communication). 

In the present paper, we further discuss the SZ ef- 
fect and address what astrophysical properties can be 
deduced with a measurement of the SZ power spectrum. 
For this, we require detailed knowledge of the covariance 
of the SZ power spectrum beyond the simple Gaussian 
sample variance. Given that the SZ effect probes the 
projected pressure distribution in the local universe, its 
statistical properties reflect those of pressure. As dis- 
cussed in detail in Cooray , the statistics of large scale 
structure pressure is highly non-Gaussian due to the asso- 
ciated non-linearities. The same non-Gaussianities lead 
to a four-point correlation function in pressure, which in 
turn, contributes non-negligibly to the power spectrum 
covariance of the SZ effect. 

In order to calculate the covariance associated with 
SZ power spectrum measurements, we extend the semi- 
analytical model presented in Cooray |^ and calculate 
the pressure trispectrum. The full covariance matrix 
now allows us to quantify the astrophysical abilities of 
SZ measurements as a probe of gas and its temperature 
properties. Previous to this study, we were unable to 
perform a similar calculation on how well the SZ effect 
probes gas and temperature properties due to the un- 
known covariance associated with the effect. Using a 
Fisher matrix approach, we quantify to what extent the 
SZ power spectrum measurements allow determination of 
various properties associated with gas and temperature 
evolution. 

Extending our calculation on the contribution of large 
scale structure gas distribution to CMB anisotropics 
through SZ effect, we also study an associated effect in- 
volving baryons associated with halos in the large scale 
structure. It is well known that the peculiar velocity of 
galaxy clusters, along the line of sight, also lead to a 
contribution to temperature anisotropics. This effect is 
commonly known as the kinetic, or kinematic, Sunyaev- 
Zel'dovich effect and arises from the baryon density mod- 
ulation of the Doppler effect associated with the velocity 
field Given that both density and velocity fields are 
involved, the kinetic SZ effect is essentially second or- 
der in density fluctuations; the thermal SZ effect is also 
second order beacuse of density and temperature depen- 
dence. Though the kinetic SZ effect was first described 
in 1^ using massive galaxy clusters, the same effect has 
been introduced under a different context by Ostriker 
and Vishniac (OV; The kinetic SZ effect can be 

considered as the OV effect extended to the non-linear 
regime of baryon fluctuations, however, it should be un- 



derstood that the basic physical mechanism responsible 
for the two effects is the same. For the purpose of this 
presentation, we will treat both OV effect and the SZ ki- 
netic effect as one contribution, though it may be easier 
to think of OV as the linear contribution while kinetic 
SZ, extending to non- linear regime will contain the total 
contribution. Such a description has been provided by 
Hu in H. 

We calculate the kinetic SZ/OV effect, hereafter sim- 
ply referred to as the kinetic SZ effect, using the model 
we developed to study the thermal SZ effect. We fur- 
ther extend this calculation to consider the correlation 
between SZ thermal and SZ kinetic effects. Since there is 
no first order cross-correlation, the lowest order contribu- 
tion to the correlation comes from a three-point function, 
or a bispectrum. Cooray & Hu discusses this bispec- 
trum in (2|. Here, we consider an additional possibility 
to measure the SZ thermal-kinetic cross-correlation via a 
two-point correlation function which involves squares of 
the temperature instead of the usual temperature itself. 
The power spectrum of squared temperatures probes one 
aspect of the trispectrum resulting from the pressure- 
baryon cross-correlation. Here, we show that there is ad- 
equate signal-to-noise for a reliable measurement of the 
SZ thermal-SZ kinetic squared power spectrum measure- 
ment in upcoming small angular scale experiments. 

The layout of the paper is as follows. In § |l| we review 
the background material relevant for current calculations 
on SZ thermal and SZ kinetic effects in the context of the 
adiabatic cold dark matter (CDM) models and using clus- 
tered dark matter halos as a way to calculate non-linear 
clustering properties. In § 3, we discuss the pressure 
trispectrum and the covariance of pressure fiuctuations 
in the local universe. In § 4, we detail the calculation 
of SZ effect due to large scale structure gas distribution 
and calculate its trispectrum. In the same section, using 
a Fisher matrix formulation, we use the full covariance, 
including non-Gaussianities, of binned SZ power spec- 
trum measurements to establish the astrophysical uses 
of the SZ effect. The kinetic SZ effect is discussed in § 5, 
and we study the cross-correlation between SZ thermal 
and SZ kinetic effects in the form of a bispectrum and 
the power spectrum of squared temperatures in § 6. We 
conclude in § 7 with a summary. 

II. GENERAL DERIVATION 

We first review the properties of adiabatic CDM mod- 
els relevant to the present calculations. We then dis- 
cuss the general properties of the halo model as applied 
to the calculation of the non-linear dark matter, baryon 
and pressure density field power spectra of the local large 
scale structure. 
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A. Adiabatic CDM Model 

The expansion rate for adiabatic CDM cosmological 
models with a cosmological constant is 

^ Hi + zf + ilKil + zf + r!A] , (1) 

where iJo can be written as the inverse Hubble distance 
today H^^ = 2997.9/i^iMpc. We follow the conventions 
that in units of the critical density 3iJo/87rG, the con- 
tribution of each component is denoted 51^, i = c for the 
CDM, g for the baryons, A for the cosmological constant. 
We also define the auxiliary quantities = r^c + ^^g and 
= 1 ~ 'Ylii which represent the matter density and 
the contribution of spatial curvature to the expansion 
rate respectively. 

Convenient measures of distance and time include the 
conformal distance (or lookback time) from the observer 
at redshift z = 



r(z) = 



dz' 



(2) 



and the analogous angular diameter distance 



1/2, 



(3) 



0, dA 



and we define r{z ~ oo) 



Note that as 
ro- 

The adiabatic CDM model possesses a two, three and 
four-point correlations of the dark matter density field as 
defined in the usual way 

(^(ki),5(k2)) = (2^)3fc(ki2)P(fcl) , (4) 
(<5(ki)<5(k2)<5(k3)) = (2^)3fc(ki23)B(ki, k2, kg) , (5) 
((5(ki)...5(k4)), = (2^f&(ki234)r(ki,k2,k3,k4), (6) 

where k^ j = k^ + . . . + kj and (5d is the delta function 
not to be confused with the density perturbation. Note 
that the subscript c denotes the connected piece, i.e. the 
trispectrum is defined to be identically zero for a Gaus- 
sian field. Here and throughout, we occasionally suppress 
the redshift dependence where no confusion will arise. 
In linear perturbation theoryfl 



2^ 



\z=0 



n+3 



(7) 



We use the fitting formulae of Eisenstein & Hu in ||2^] 
to evaluate the transfer function T{k) for CDM models. 



Here, Sh is the amplitude of present-day density fluctu- 
ations at the Hubble scale. 

The bispectrum in perturbation theory is given by ^ 

B'"\kp, k„ k,) = 2F|(kp, k,)Pikp)P{k,) + 2 Perm. , 

(8) 

with F2 term given by second order gravitational pertur- 
bation calculations. 

Similarly, the perturbation theory trispectrum follows 
from Fry p6[ : 

T"" = 6[^^3^(ki,k2,k3)F(fci)P(fc2)F(fc3) +Perm.] 
+ 4 [F|(ki2, -ki)F|(ki2, k3)P(fci)P(fci2)F(fc3) + Perm.] . (9) 

The permutations involve a total of 4 terms in the first 
set and 12 terms in the second set. 

In linear theory, the density field may be scaled back- 
wards to higher redshift by the use of the growth function 
G{z), where 6{k,r) = G(r)(5(fc,0) 



G(r) oc 



iJ(r) 



Jz{r) 



dz'{l + z') 



Hn 



H{z') 



(10) 



Note that in the matter dominated epoch G cx a = (1 + 

The continuity equation relates the density and veloc- 
ity fields, in the linear regime, via. 



'iG5{k,Q) 



k^ 



(11) 



where overdots represent derivatives with respect to ra- 
dial distance r. In this paper, we ignore the contribution 
to the velocity field within individual non-linear struc- 
tures, including a curl component, and only consider this 
potential bulk flow. This large scale velocity held is inde- 
pendent of the non-linear density fluctuations, and thus, 
there is no correlation between this bulk flow and the 
density field within halos. 

For fluctuation spectra and growth rates of interest 
here, reionization of the universe is expected to occur 
rather late Zri ^ 50 such that the reionized media is opti- 
cally thin to Thomson scattering of CMB photons t ^ 1. 
The probability of last scattering within dr of r (the vis- 
ibility function) is 



g^ie-^ = XHoTHil + z)h- 



Here r(r) = /J drr is the optical depth out to r, X is 
the ionization fraction, 



•'•It should be understood that "lin" denotes here the lowest 
non-vanishing order of perturbation theory for the object in 
question. For the power spectrum, this is linear perturbation 
theory; for the bispectrum, this is second order perturbation 
theory, etc. 



^Tho kernels are derived in |25| (see, equations A2 and 
A3 of |2^]; note that their Pn = Fn), and we have written 
such that the symmetric form of FnS are used. The use of 
the symmetric form accounts for the factor of 2 in Eqs. (^) 
and factors of 4 and 6 in (19^. 
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FIG. 1. The dark matter (NFW) profile and the ones pre- 
dicted by the hydrostatic equihbrium for g ctS , ctS 3i function 
of the 6 parameter (see, Eq. M) with = 100. For compari- 
son, we also show a typical example of the so-called (3 model 
(1 + /r^)~^^^'^ which is generally used as a fitting function 
for X-ray and SZ observations of clusters. 



TH = 0.0691(1 - Yp 



(13) 



is the optical depth to Thomson scattering to the Hubble 
distance today, assuming full hydrogen ionization with 
primordial helium fraction of 1^. Note that the ionization 
fraction can exceed unity: X = {1 — 3Yp/4)/{l — Yp) for 
singly ionized helium, X = {1 ~ Yp/2) /{I — Yp) for fully 
ionized helium. 

Although we maintain generality in all derivations, we 
illustrate our results with the currently favored ACDM 
cosmological model. The parameters for this model are 
= 0.30, fig = 0.05, nA = 0.65, h = 0.65, Yp = 0.24, 
n = 1, X = 1, with a normalization such that mass fluc- 
tuations on the 8h Mpc~^ scale is erg = 0.9, consistent 
with observations on the abundance of galaxy clusters 
psf and COBE normalization |2^. A reasonable value 
here is important since higher order correlations is non- 
linearly dependent on the amplitude of the density field. 

For the Sunyaev-Zel'dovich effects discussed here, we 
are more interested in the clustering properties of pres- 
sure and baryons, rather than the dark matter density 
field. We do not have a reliable way to calculate the 
pressure, baryon power spectra and their higher order 
correlations and cross-correlations analytically. We will 
introduce the semi-analytic halo model for this purpose 
following pi]. 



B. Halo Approach 

In order to calculate the non-linear clustering of the 
dark matter, baryon and pressure density fields, we use 



the halo model described in ||30|-p2[. 

Underlying the halo approach is the assertion that dark 
matter halos of virial mass M are locally biased tracers of 
density perturbations in the linear regime. In this case, 
functional relationship between the over-density of halos 
and mass can be expanded in a Taylor series 



4(x,M;z) = 

bo + bi{M;z)S ) + -b2{M; z)6fi^{x; z) + . . . (14) 



The over-density of halos can be related to more famil- 
iar mass function and the halo density profile by noting 
that we can model the fully non-linear density field as a 
set of correlated discrete objects or halos with profiles 



P(x) 



Ph(x - Xj; Af^; 



(15) 



where the sum is over all halo positions. The density 
fluctuation in Fourier space is 



4(k,M,). 



(16) 



Following Peebles in |g7|, we divide space into suffi- 
ciently small volumes 6V that they contain only one or 
zero halos of a given mass and convert the sum over halos 
to a sum over the volume elements and masses 

S{k)= J2 me*-"^(5h(k,Mi). (17) 

Vi,A/i 

By virtue of the small volume element rii = nf ~ n'^ = 
1 or following 

Using the fact that halos are biased tracers of the den- 
sity field such that their number density fluctuates as 



dn 
dM 



(x;z) 



dn 
dM 



{M; z) 



X [bo + 6i(A/; z)Jun(x; z) + -b2{M; z)Sl^{x; z) 



we can write 



(18) 



(19) 



dMi 

(nms) = (ni)^ 5i2 + (ni) (na) [6^ + 6i(Afi)6i(M2) , 

X ((5lin(xi)(5iin(x2))] . 

{nin2n-i) = . . . . (20) 

In Eq. |l^, 6o = 1 and 5i2 is the Dirac delta function and 
we have included only the lowest order terms. The halo 
bias parameters are given in 133]: 



6i(Af;z) = 1 + 



iy^{M;z) ~ I 



(21) 



and 
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&2(Af;z) = — [6i(M;z)-l] + 



iy^{M;z) - 3 
(7^{M;z) 



(22) 



Here, j/(M; z) = 5c/(t{M; z) and cr(M; z) is the rms fluc- 
tuation within a top-hat filter at the virial radius corre- 
sponding to mass M, and 6c is the threshold over-density 
of spherical collapse (see |Q for useful fitting functions) . 

The derivation of the higher point functions in Fourier 
space is now a straightforward but tedious exercise in 
algebra (see, ]49|]). The Fourier transforms inherent in 
eqn. ( p^ ) convert the correlation functions in eqn. ( pO| ) 
into the power spectrum, bispectrum, trispectrum, etc., 
of perturbation theory. 

Following and [Q, it is now convenient to define 
a general integral over the halo mass function dn/dM. 
Though we presented the description of halo clustering 
for dark matter, we can generalize this discussion to con- 
sider any physical property associated with halo; one sim- 
ply relates the over-density of halos through the density 
profile corresponding to the property of interest in Eq. nj. 
Since we will encounter dark matter, pressure and baryon 
density fields through out this paper, we write a general 
integral that applies to all these three properties 



diT 

{ki,...,k^;z)EE I dM—{M,z)bp{M;z) 



xT,{M;z)^y,,{h,M;z) 



dM 

y,^{k^,M;z) . 



(23) 



Here, in addition to the dark matter, to account for clus- 
tering properties of pressure associated with baryons in 
large scale structure, we have introduced the electron 
temperature, Tc{M]z). 

In Eq. |2^, the three-dimensional Fourier transform of 
the density fluctuation, through the halo profile of the 
density distribution of any physical property Pi(r, M; z), 
is 



yi{k,M;z) = — 
Pbi 



dr 4:Trr^pi{r, M; z) 



sin(fcr) 
kr 



(24) 



with the background mean density of the same quantity 
given by pbi. Since in this paper we discuss the dark 
matter, pressure and baryons, the index i will be used 
to represent the density, 6 (with y = ys), the baryons, g 
(with y = yg), or pressure, H (with = yg). Note that for 
both dark matter density and baryon clustering, -q — Q, 
as there is no temperature contribution, but for cluster- 
ing of pressure, rj = p when ii . . .i^ describes pressure; 
therefore, we will no longer make use of the index 77 in 
the following discussion. One additional note here is that 
the profile used for baryons will be the same as the pro- 
file that we will use for pressure. The only difference be- 
tween baryon clustering and pressure clustering is that we 
weigh the latter with the electron temperature, leading 
to a selective contribution from electrons with the highest 
temperature, while the former includes all baryons. 
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FIG. 2. The pressure power spectrum (a) and 
square-configuration trispectrum (b) today (z — 0) broken 
into individual contributions under the halo description. The 
line labeled 'bias' shows the pressure bias relative to the dark 
matter power spectrum under the halo model. In (b), we 
show the square configuration trispectrum (see text). In both 
(a) and (b), at all scales relevant for the Sunyaev-Zel'dovich 
effect, the single halo term dominates. 



C. Ingredients 

The dark matter profile of collapsed halos are taken to 
be the NFW pS] with a density distribution 



Ps{r) 



{r/rs)il + r/r,y 



(25) 



The density profile can be integrated and related to the 
total dark matter mass of the halo within r„ 



Ms = Anpsr^ 



log(l +c) 



1 



(26) 



where the concentration, c, is ry/rg. Choosing as the 
virial radius of the halo, spherical collapse tells us that 
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FIG. 3. The mass dependence on the dark matter power 
spectrum (a) and pressure power spectrum (b). Here, we 
show the total contribution broken in mass limits as written 
on the figure. As shown in (a), the large scale contribution to 
the dark matter power comes from massive halos while small 
mass halos contribute at small scales. For the pressure, in 
(b), only massive halos above a mass of 10^ Mq contribute 
to the power. 



M = 47rr^A(z)pb/3, where A(z) is the over-density of 
collapse and pt is the background matter density today. 
We use comoving coordinates throughout. By equating 
these two expressions, one can eliminate and describe 
the halo by its mass M and concentration c. We use 
the Press-Schechter mass function to describe the 
halo mass function written as dn/dM{M\z). Following 
the results from ACDM simulations by [g^j we take a 
concentration-mass relationship which is consistent with 
the mean: 



c(M; z) = 9(1 -hz)"^ 



M 



n -0.13 



M4z) 



(27) 



where Mi,{z) is the non-linear mass scale at which the 
peak-height threshold v{M] z) is equal to one. Due to 
computational limitations, we ignore the distribution of 
concentrations, for a given mass, observed in simulations 
by and only use the mean value. As discussed in [ pO[ , 
the concentration distribution is important for higher 
moments and ignoring the distribution in present work, 
it is likely that we underestimate the non-Gaussian con- 
tribution from the single halo term at a level of 20%. We 
also expect uncertainties at the same level due to our 
modeling of halos as smooth spherical distributions and 
by ignoring the substructure and asphericity of halos. 

The gas density profile, Pg{r), is calculated assuming 
hydrostatic equilibrium between the gas distribution and 
the dark matter density field within a halo. This is a 
valid assumption given that current observations of ha- 
los, mainly galaxy clusters, suggest the existence of regu- 
larity relations, such as size-temperature (e.g., [Q), be- 
tween physical properties of dark matter and baryon dis- 
tributions. Using the hydrostatic equilibrium allows us 
to write the baryon density profile, Pg{r), within halos as 



1 



hvs j T 



where 6 is a constant, for a given mass, 
^ _ ATTGfj.mpPsr^ 



(28) 



(29) 



with the Boltzmann constant, [£9j and p = 0.59, cor- 
responding to a hydrogen mass fraction of 76%. Equa- 
tion |2^ is derived only under the assumption of hydro- 
static equilibrium for the gas distribution in a dark mat- 
ter profile given by the NEW equation. The hydrostatic 
equilibrium allows a basic physical assumption to relate 
pressure and dark matter. Though effects such as merg- 
ers between halos is expected to result in a violation of 
this assumption, we do not consider variants since pre- 
vious predictions made with the above halo have been 
found to be consistent with numerical simulations by Re- 
fregier & Teyssier ]l7t . In above, the normalization pgo is 
determined under the assumption of a constant gas mass 
fraction for halos comparable with the universal baryon 
to dark matter ratio: fg = Mg/Mg = Vlg/Q,m- When 
investigating astrophysical uses of the SZ effect, we will 
vary this parameter and consider variations of gas frac- 
tion as a function of mass and redshift. 

The electron temperature can be calculated based on 
the virial theorem or similar arguments as discussed in 
Ip^. Using the virial theorem, we can write 



jGpnipMs 
3r„ 



(30) 



with 7 = 3/2. Since cx Mg^'^{l + z) ^ in physical coor- 
dinates, Te cx M^/^{l + z). The average density weighted 
temperature is 



6 



{T^), = J dM y^^{M, z)T,{M, z) . 



(31) 



The total gas mass present in a dark matter halo within 
r,, is 



(32) 



In Fig. |l|, we show the NFW profile for the dark mat- 
ter and arbitrarily normalized gas profiles predicted by 
the hydrostatic equilibrium and virial theorem for several 
values of b. As b is decreased, such that the temperature 
is increased, since b cx 1/Te, the turn over radius of the 
gas distribution shifts to higher radii. The relative nor- 
malization of the gas profiles are set with a gas fraction 
value of 0.1, while the NFW profile is arbitrarily nor- 
malized with ps — 1; the gas profiles, therefore, simply 
scale with the same factor related to true value of ps ■ As 
an example, we also show the so-called /3 model that is 
commonly used to describe X-ray and SZ observations of 
galaxy clusters and for the derivation purpose of the Hub- 
ble constant by combined SZ/X-ray data. The (3 model 
describes quite accurately the underlying gas distribution 
predicted by the gas profile in Eq. |2^, though, we find dif- 
ferences at the outermost radii of halos. This difference 
can be used as a way to establish the hydrostatic equilib- 
rium for clusters with respect to a NFW like dark matter 
distribution, though, any difi^erence of gas distribution at 
large radii may need to be accounted in the context of 
possible substructure and mergers. 

A discussion on the comparison between the gas profile 
used here and the /3 model is available in In addi- 
tion, we refer the reader to [Q for a full detailed dis- 
cussion on issues related to modeling of pressure power 
spectrum using halos and associated systematic errors. 
Comparisons of the halo model predictions with numer- 
ical simulations are available in and ||l^. Similarly, 
issues related to modeling of the dark matter clustering 
using halos is discussed in for the bispectrum and 
pO| for the trispectrum. 



D. Power Spectra and Trispectra 

For the calculations presented in this paper, we will 
encounter the power spectrum and trispectrum involv- 
ing dark matter, pressure and baryons and their cross- 
correlations. We now write down the relevant Fourier 
space correlations under the halo approach. To general- 
ize the discussion, we will use the index i to represent the 
property of interest. 

In general, the power spectrum of these three quanti- 
ties under the halo model now becomes iSOtl: 



Pl\k)^ll,{k,k), 
Pl\k)^[ll{kfp'^-{k), 



(33) 
(34) 

(35) 



where the two terms represent contributions from two 
points in a single halo (Ih) and points in two different 
halos (2h) respectively. 

Similar to above, we can also define the cross power 
spectra between two fields i and j as 



-.,ik)^P^{k)+P^\k). 



i^f(fc) = /i,(fc)/i,.(fc)F"" 



(fc). 



(36) 
(37) 
(38) 



It is also useful to define the bias of one field relative to 
the dark matter density field as 



biasi(fc) = 



Ps{k) 



(39) 



We can also define a dimensionless correlation coefficient 
between the two fields as 



P,,{k) 



(40) 



During the course of this paper, we will encounter, and 
use, cross-power spectra such as the one involving baryon 
and pressure, Pgn, and dark matter and pressure. Pan- 

As described in |2^, the contributions to the trispec- 
trum may be separated into those involving one to four 
halos 



T,, = T. 



Ih 



2h 



T, 



3h 



Ah 



(41) 



where here and below the argument of the trispectrum is 
understood to be (ki, k2, ka, k4). The term involving a 
single halo probes correlations of the property of interest, 
e.g., dark matter, within that halo 



T 



Ih 



j(fcl,fc2,A:3,^4) 



(42) 



and is independent of configuration due to the assumed 
spherical symmetry for our halos. 

The term involving two halos can be further broken up 
into two parts 



rp2h rp2h 



T. 



2h 



(43) 



which represent taking three or two points in the first 
halo 



ah 



^31 



rp2h pi: 



%k,)il,,,{k2MM)iUk,)- 
\k,2)ilAki,k2)iiu{hM) 



3 Perm., (44) 
-2 Perm. (45) 



The permutations involve the 3 other choices of ki for the 
ll j term in the first equation and the two other pairings 
of the kiS for the l\ terms in the second. Here, we have 
defined ki2 = ki -I- k2; note that ki2 is the length of one 
of the diagonals in the configuration. 

The term containing three halos can only arise with 
two points in one halo and one in each of the others 
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T^'^ = i?""(ki,k2,k34)4,,(fc3,fc4)/i^,,(fcl)/L(fc2) 

+P'"^h)P'"^k2)llu{h, h)ll{h)ll{k2) + 5 Perm. , (46) 

where the permutations represent the unique pairings of 
the fci's in the l2,ii factors. This term also depends on 
the configuration. 

Finally for four halos, the contribution is 



T 



Ah 



I 



U''i)iUk2)iUh)ii{h){ 



rpl] 



xP""(fci)P""(fc2)P""(fc3) + 3 Perm.] }, (47) 

where the permutations represent the choice of ki in the 
in the brackets. Because of the closure condition 
expressed by the delta function, the trispectrum may be 
viewed as a four-sided figure with sides k.^. It can alter- 
nately be described by the length of the four sides ki plus 
the diagonals. We occasionally refer to elements of the 
trispectrum that differ by the length of the diagonals as 
different configurations of the trispectrum. In the rest 
of the paper, we will encounter the pressure trispectrum 
Tn = Tnnnn and the pressure-baryon cross trispectrum: 



III. PRESSURE TRISPECTRUM AND POWER 
SPECTRUM COVARIANCE 

Following we can relate the trispectrum to the 
variance of the estimator of the binned power spectrum 



1 

V 



(Pk 



<5^,(-k)<5n(k), 



(48) 



where the integral is over a shell in fc-space centered 
around /ci, Vsi ~ AirkfSk is the volume of the shell 
and V is the volume of the survey. Recalling that 
(5(0) V/{2Tr)^ for a finite volume. 



pp. 



P ) P. 



1 

V 



(2^)^ 



where 



Th 



2PfS, 



Tn(k, 



(49) 



-k,) . (50) 



Notice that though both terms scale in the same way with 
the volume of the survey, only the Gaussian piece nec- 
essarily decreases with the volume of the shell. For the 
Gaussian piece, the sampling error reduces to a simple 
root-N mode counting of independent modes in a shell. 
The trispectrum quantifies the non-independence of the 
modes both within a shell and between shells. Calcu- 
lating the covariance matrix of the power spectrum esti- 
mates reduces to averaging the elements of the trispec- 
trum across configurations in the shell. 



A. Discussion 

In Fig. |^(a), we show the logarithmic power spec- 
trum of pressure and dark matter such that A^(fc) = 
k^P{k)/2TT^ with contributions broken down to the Ih 
and 2h terms today. As shown, the pressure power spec- 
trum depicts an increase in power relative to the dark 
matter at scales out to few h Mpc~^, and a decrease 
thereafter. The decrease in power at small scales can 
be understood through the relative contribution to pres- 
sure as a function of the halo mass. In Fig. ||, we break 
the total dark matter power spectrum (a) and the to- 
tal pressure power spectrum (b), to a function of mass. 
As shown, contributions to dark matter come from mas- 
sive halos at large scales and by small mass halos at 
small scales. The pressure power spectrum is such that 
through temperature weighing, with Tf, oc M^/^ depen- 
dence, the contribution from low mass halos to pressure 
is suppressed relative to that from the high mass end. 
Thus, the pressure power spectrum, at all scales of in- 
terest, can be described with halos of mass greater than 
10^^ Mq. a comparison of the dark matter and pres- 
sure power spectra, as a function of mass, in Fig. |^ re- 
veals that the turn-over in the pressure power spectrum 
results at an effective scale radius for halos with mass 
greater than 10'^'' Mq. As we find later, this effective 
turn-over produces a distinct signature in the SZ angular 
power spectrum inolving a turnover of SZ power at mul- 
tipoles of few thousand; The presence of this turnover 
and the exact multipole range can be used as a probe of 
gas and temperature physics. 

For the covariance of pressure, and also for the covari- 
ance of the SZ power spectrum, we are mainly interested 
in terms of the pressure trispectrum involving configura- 
tions that result in Tn(ki, — ki, k2, — k2), i.e. parallelo- 
grams which are defined by either the length fci2 or the 
angle between ki and k2. These are the configurations 
that contribute to the power spectrum covariance. For 
illustration purposes we will take fci = fc2 and the angle 
to be 90° (k2 = kj_) such that the parallelogram is a 
square. It is then convenient to define 



ALq(fc)^ |^<'(k,-k,k^,-kx), 



(51) 



such that this quantity scales roughly as the logarithmic 
power spectrum itself A^(fc). This spectrum is shown in 
Fig. ^(b) with the individual contributions from the Ih, 
2h, 3h, 4h terms shown. As shown, almost all contribu- 
tions to the pressure trispectrum come from the single 
halo term. 

Using the pressure trispectrum, we can now predict 
the pressure covariance and, more appropriately, corre- 
lations in the binned measurements of the pressure. The 
predictions made here with the halo model to describe 
pressure can easily be tested in numerical simulations 
and the accuracy of the halo model can be further stud- 
ied. For this purpose, we calculate the covariance matrix 
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TABLE I. 

Pressure Power Spectrum Correlations 



k 


0.031 


0.044 


0.058 


0.074 


0.093 


0.110 


0.138 


0.169 


0.206 


0.254 


0.313 


0.385 


0.031 


1.000 


0.243 


0.337 


0.386 


0.398 


0.396 


0.385 


0.371 


0.359 


0.350 


0.343 


0.337 


0.044 


(0.019) 


1.000 


0.442 


0.546 


0.576 


0.578 


0.566 


0.547 


0.533 


0.521 


0.513 


0.507 


0.058 


(0.041) 


(0.036) 


1.000 


0.653 


0.728 


0.745 


0.736 


0.715 


0.699 


0.687 


0.678 


0.672 


0.074 


(0.065) 


(0.075) 


(0.062) 


1.000 


0.807 


0.865 


0.868 


0.847 


0.832 


0.819 


0.811 


0.805 


0.093 


(0.086) 


(0.111) 


(0.118) 


(0.102) 


1.000 


0.898 


0.932 


0.918 


0.902 


0.890 


0.882 


0.876 


0.110 


(0.113) 


(0.153) 


(0.183) 


(0.189) 


(0.160) 


1.000 


0.955 


0.959 


0.946 


0.934 


0.926 


0.921 


0.138 


(0.149) 


(0.204) 


(0.255) 


(0.299) 


(0.295) 


(0.277) 


1.000 


0.978 


0.973 


0.961 


0.953 


0.947 


o.iby 


(0.1 1 1) 




/A QAO\ 

(O.oUzj 


/"A OCiQ\ 

(O.oboj 


/A oc:o\ 

(U.oboj 


/'A AOO\ 

(0.4o3j 


/A AO A\ 

(U.4J4J 


1 AAA 
i.UUU 


A AOO 

0.982 


A ATO 

0.9(2 


A ACO 

0.962 


A A K^' 


0.206 


(0.186) 


(0.261) 


(0.334) 


(0.334) 


(0.412) 


(0.541) 


(0.580) 


(0.592) 


1.000 


0.985 


0.973 


0.965 


0.254 


(0.186) 


(0.264) 


(0.341) 


(0.341) 


(0.425) 


(0.576) 


(0.693) 


(0.737) 


(0.748) 


1.000 


0.986 


0.973 


0.313 


(0.172) 


(0.251) 


(0.328) 


(0.328) 


(0.412) 


(0.570) 


(0.698) 


(0.778) 


(0.839) 


(0.858) 


1.000 


0.986 


0.385 


(0.155) 


(0.230) 


(0.305) 


(0.305) 


(0.389) 


(0.549) 


(0.680) 


(0.766) 


(0.848) 


(0.896) 


(0.914) 


1.000 


A^ 


1.04 


1.13 


1.29 


1.64 


1.95 


2.24 


3.58 


6.09 


9.27 


16.4 


21.2 


29.8 




1.00 


1.01 


1.02 


1.03 


1.04 


1.07 


1.14 


1.23 


1.38 


1.61 


1.90 


2.26 



NOTES. — Diagonal normalized covariance matrix of the binned pressure (upper triangle) and dark matter density field (lower triangle) 
pow er spectrum with k values in units of h Mpc^^. Lower triangle (parenthetical numbers) displays the covariance from halo model in 
|20[ . Final line shows the fractional increase in the errors (root diagonal covariance) due to non-Gaussianity as calculated under the halo 
model for pressure and dark matter density field. 



Cy from Eqn. (pQ) with the bins centered at ki and vol- 
ume Vsi — 'inkfSki. The binning scheme used here is 
the one we utihzed in |^ to calculate the binned dark 
matter power spectrum correlations. 

In Table |, we tabulate the pressure (upper triangle), 
and for comparison dark matter (lower triangle), corre- 
lation coefRcients 

a, = -^i= . (52) 

The dark matter correlations are from the halo based 
predictions by There, for the dark matter, we sug- 
gested that the halo model predicted correlations agree 
with numerical simulations of jiot typically better than 
±0.1, even in the region where non-Gaussian effects dom- 
inate, and that the the qualitative features such as the 
increase in correlations across the non-linear scale are 
preserved. As we do not have measurements of the pres- 
sure correlations from simulations, we cannot perform a 
detailed comparison on the accuracy of the halo model 
predictions for pressure here. 

A further test on the accuracy of the halo approach 
is to consider higher order real-space moments such as 
the skewness and kurtosis. In [|l8|, we discussed the SZ 
skewness under the halo model. As discussed in detail by 
halo model predictions agree remarkably well with 
numerical simulations, especially for the pressure and SZ 
power spectra, though, detailed comparisons still remain 
to be made with respect to bispectrum and trispectrum. 

Even though the dark matter halo formalism provides 
a physically motivated means of calculating the statistics 
of the dark matter density field and associated properties 



such as pressure, there are several limitations of the ap- 
proach that should be borne in mind when interpreting 
results. The approach assumes all halos share a parame- 
terized spherically-symmetric profile and this assumption 
is likely to affect detailed results on the configuration de- 
pendence of the bispectrum and trispectrum. Since we 
are considering a weighted average of configurations, our 
predictions presented here may be insufficient to estab- 
lish the validity of the trispectrum modeling in general. 
Further numerical work is required to quantify to what 
extent the present approach reproduces simulation re- 
sults for the full trispectrum. We do not consider such 
comparisons here, other than to suggest that the halo 
model has provided, at least qualitatively, a consistent 
description better than any of the arguments involving a 
biased description of gas tracing the dark matter etc. 

IV. SZ THERMAL EFFECT 

The Sunyaev-Zel'dovich (SZ; §]) effect arises from the 
inverse- Compton scattering of CMB photons by hot elec- 
trons along the line of sight. The temperature decrement 
along the line of sight due to SZ effect can be written as 
the integral of pressure along the same line of sight 

y ^ 7^ = g{x) f dra{r)'^n,{r)T,{r) (53) 

J CMB J rUeC-' 

where ax is the Thomson cross-section, Ug is the elec- 
tron number density, r is the comoving distance, and 
g{x) = a;coth(a;/2) — 4 with x = hv/ksTcMB is the spec- 
tral shape of SZ effect. At Rayleigh-Jeans (RJ) part of 
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10'' 



10'" 
10" 

gio^^ 

10^' 
10" 
10"' 

10' 10' 10" 10' 10" 10' 

/ / 

FIG. 4. The angular power spectra of SZ thermal and kinetic effects. As shown in (a), the thermal SZ effect is dominated by 
individual halos, and thus, by the single halo term, while the kinetic effect is dominated by the large scale structure correlations 
depicted by the 2-halo term. In (b), we show the mass dependence of the SZ thermal and kinetic effects with a maximum mass 
of 10^® and 10^'^ ^Q- The SZ thermal effect is strongly dependent on the maximum mass, while due to large scale correlations, 
kinetic effect is not. 




the CMB, g{x) = —2. For the rest of this paper, we 
assunie observations in the Rayleigh- Jeans regime of the 
spectrum; an experiment such as Planck with sensitivity 
beyond the peak of the spectrum can separate out SZ 
contributions based on the spectral signature, g{x) [|2|. 



A. Power Spectrum and Trispectrum 

The SZ power spectrum and trispectrum are defined 
in the flat sky approximation in the usual way 

(y(li)y(l2)) =(2^)^fc(li2)CP, 
(y(li) . . . y(l4)), = (2^)2fc(li234)rSZ(li, l2, 13, 14) . (54) 

These can be written as a redshift projection of the pres- 
sure power spectrum 



SZ 



T^^ = I dr 



''A 



dA' 

il il il 



(55) 
(56) 



where dA is the angular diameter distance. At RJ part 
of the frequency spectrum, the SZ weight function is 



(57) 



where is the mean electron density today. In deriving 
Eq. (|5^) , we have used the Limber approximation |Q by 
setting k — l/dA and flat-sky approximation. 

In Fig. ^(a), we show the SZ power spectrum due to 
baryons present in virializcd halos. As shown, most of 
the contributions to SZ power spectrum comes from in- 
dividual massive halos, while the halo-halo correlations 
only contribute at a level of 10% at large angular scales. 
This is contrary to, say, the lensing convergence power 
spectrum discussed in Q, where most of the power at 
large angular scales is due to the halo-halo correlations. 
The difference can be understood by noting that the SZ 
effect is strongly sensitive to the most massive halos due 
to T oc M^/^ dependence in temperature and to a lesser, 
but somewhat related, extent that the weight function 
increases towards low redshifts. The lensing weight func- 
tion selectively probes the large scale dark matter den- 
sity power spectrum at comoving distances half to that 
of background sources (z ~ 0.2 to 0.5 when sources are 
at a redshift of 1), but has no extra dependence on mass. 
The fact that the SZ power spectrum results mainly from 
the single halo term also results in a sharp reduction of 
power when the maximum mass used in the calculation is 
varied. For example, as discussed in ||l^ and illustrated 
in Fig 0(b), with the maximum mass decreased from 10^^ 
to 10^ Mq, the SZ power spectrum is reduced by a fac- 
tor nearly two orders of magnitude in large scales and an 
order of magnitude at Z ~ 10**. 

In Fig ^, we show the scaled trispectrum 



10 



10" 



10" 



^ 10-' 



10" 



1h 



2h 



3h 
4h 



10' 



10' 



/ 



FIG. 5. SZ trispectrum for a configuration defined by a 
square under the halo modeL As shown, single halo term 
dominates the contribution to trispectrum at all multipoles 
ranging from large angular scales to small angular scales. The 
dependence on the single halo term is consistent with the gen- 
eral non-Gaussian behavior of the SZ effect and it is significant 
non-linearity. 



. SZ 



(0 = (1' 



'1 



Nl/3 



(58) 



where — I and 1-1^ = 0. The projected SZ trispectrum 
again shows the same behavior as the pressure trispec- 
trum. As shown, the contributions to the trispectrum 
essentially comes from the single halo term at all mul- 
tipoles. This is consistent with our observation that SZ 
power spectrum is essentially dominated by the corre- 
lations of pressure within halos. As discussed in [^ , 
the SZ bispectrum is also dominated by the single halo 
term. Given this dependence on the single halo term, 
for the rest of the discussion involving SZ covariance, we 
will only use the single halo term and ignore the con- 
tributions arising from large scale correlations associated 
with halos. 



B. SZ Power Spectrum Covariance 

For the purpose of this calculation, we assume that 
planned small angular scale SZ experiments will measure 
the thermal SZ power spectrum in bins of logarithmic 
band powers at several U 's in multipole space with thick- 
ness 5li. Thus, we can write the band power measure- 
ments as 



C^ 



si ^si 27r 



y(i)y(-i): 



(59) 



where As{li) = J cPl is the area of 2D shell in multipole 
and can be written as As{li) = 2Trli6li + Tr{Sli)'^ . 
We can now write the signal covariance matrix as 



10"' 



10 ■ 



10 



Gaussian 

with non-Gaussianities 



10' 



10^ 



10' 



/ 



FIG. 6. The fractional errors in the measurements of the 
SZ band powers. Here, we show the fractional errors un- 
der the Gaussian approximation, and the total including 
non-Gaussianities. As shown, the total contribution as a func- 
tion of mass is sensitive to the presence of most massive halos 
in the universe. The non-Gaussian term is essentially domi- 
nated by the single halo term. 



2Cf + T; 



SZ 



■SZ 



A-si 



2;2 



T^^(]- —1 1 —1 



(60) 



where A is the area of the survey in steradians. Again 
the first term is the Gaussian contribution to the sample 
variance and the second term is the non-Gaussian con- 
tribution. A realistic survey will also have an additional 
noise variance due to the instrumental effects and a co- 
variance resulting from the uncertainties associated with 
the separation of the SZ effect from thermal CMB and 
other foregrounds. 

In Fig. ^, we show the fractional error. 



(61) 



for bands h given in Table H following the binning scheme 
used in pol] for the weak lensing power spectrum. 

In Fig. Id, the dashed line shows the Gaussian error 
while the solid line shows the total covariance with the 
addition of the SZ trispectrum (Eq. ^0|). At all multi- 
poles, the non-Gaussianities from the trispectrum domi- 
nates the variance. As we discussed for the power spec- 
trum, however, a reduction in the maximum mass of the 
halos used for the SZ calculation leads to a sharp de- 
creases in the non-Gaussianities. With a mass cut at 
10^^ Mq, shown by the dotted line, we see that the total 
variance is consistent with the Gaussian variance out to 
I - 1000. 

We can now write the correlation between the bands 



as 



C,;. 



\/ Cii C. 



(62) 
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TABLE II. 

SZ Thermal Power Spectrum Correlations 





97 


138 


194 


271 


378 


529 


739 


1031 


1440 


2012 


2802 


3905 


5432 


7568 


97 


1.00 


0.90 


0.71 


0.49 


0.30 


0.18 


0.09 


0.05 


0.02 


0.01 


0.00 


0.00 


0.00 


0.00 


138 


(0.00) 


1.00 


0.90 


0.70 


0.49 


0.30 


0.17 


0.09 


0.05 


0.03 


0.02 


0.02 


0.01 


0.00 


194 


(0.00) 


(0.00) 


1.00 


0.89 


0.70 


0.48 


0.30 


0.17 


0.09 


0.05 


0.03 


0.03 


0.01 


0.00 


271 


(0.00) 


(0.00) 


(0.00) 


1.00 


0.89 


0.69 


0.47 


0.30 


0.17 


0.09 


0.05 


0.02 


0.01 


0.00 


378 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


1.00 


0.89 


0.69 


0.49 


0.30 


0.17 


0.09 


0.05 


0.03 


0.01 


529 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


1.00 


0.88 


0.69 


0.48 


0.30 


0.17 


0.17 


0.11 


0.07 


739 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


1.00 


0.88 


0.68 


0.49 


0.38 


0.28 


0.20 


0.13 


1031 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


1.00 


0.87 


0.68 


0.56 


0.43 


0.31 


0.21 


1440 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


1.00 


0.88 


0.69 


0.61 


0.45 


0.31 


2012 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.01) 


(0.05) 


1.00 


0.88 


0.69 


0.60 


0.42 


2802 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.02) 


(0.09) 


(0.39) 


1.00 


0.88 


0.70 


0.56 


3905 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.02) 


(0.08) 


(0.36) 


(0.84) 


1.00 


0.87 


0.70 


5432 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.01) 


(0.06) 


(0.29) 


(0.65) 


(0.86) 


1.00 


0.88 


7568 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.00) 


(0.01) 


(0.04) 


(0.20) 


(0.53) 


(0.70) 


(0.88) 


1.00 



NOTES. — Covariance of the binned power spectrum for the SZ effect. Upper triangle displays the covariance found when a perfect 
frequency cleaned SZ map is used to determine the SZ power spectrum. Lower triangle (parenthetical numbers) displays the covariance 
found when the variance is dominated by the primary anisotropy contribution, as in a measurement of the SZ power spectrum in a CMB 
primary fluctuations dominated map. 



In Table || we tabulate the SZ correlations under the 
assumption that the SZ power spectrum is measured in- 
dependently, say in a frequency cleaned map, (upper tri- 
angle) and is measured in the CMB primary dominated 
map (lower triangle). The correlations along individual 
columns increase (as one goes to large Vs or small angular 
scales) and the maximum values are reached at / 5000 
consistent with the general behavior of the trispectrum. 

In Fig. 1^, we show the non-Gaussian trispectrum cor- 
relation coefficient given by 

Cr - ■ (63) 

\/ J- 21 -1-2 j 

As shown here, the increase in non-Gaussian correlation 
is consistent with the fact that at all scales it is the single 
halo term which dominates the non-Gaussian contribu- 
tion. Since the power spectrum is dominated by correla- 
tions in single halos, the fixed profile of the halos correlate 
the power in all the modes and the correlations between 
adjacent modes are significant. 

The calculation, or experimental measurement, of the 
full SZ covariance is necessary for the interpretation of 
observational results on the power spectrum. The up- 
coming SZ surveys, where the power spectrum will be 
measured, is likely to be limited to a small area on the 
sky. Thus, in the absence of many fields where the covari- 
ance can be estimated directly from the data, the halo 
model based approach suggested here provides a useful, 
albeit model dependent, quantification of the covariance. 
As suggested for weak lensing observations in |^l| and 
discussed in [ po[ , as a practical approach one could imag- 
ine taking the variances estimated from the survey under 
a Gaussian approximation, after accounting for uneven 



sampling and edge effects, and scaling it up by the non- 
Gaussian to Gaussian variance ratio of the halo model 
along with inclusion of the band power correlations. Ad- 
ditionally, using the covariance as the one calculated here, 
one can use the approach well known in the fields of CMB 
and galaxy power spectrum measurements to decorrelate 
band powers . 

C. Astrophysical Uses of the SZ Power Spectrum 

The calculation of the full covariance matrix now al- 
lows us to study how well the SZ power spectrum mea- 
sures certain astrophysical and cosmological parame- 
ters. The upcoming CMB power spectrum measure- 
ments, complemented by the related local universe ob- 
servations such as the galaxy power spectrum or super- 
novae, are expected to constrain most of the cosmological 
parameters to a reasonable accuracy |Q. Thus, we ig- 
nore the possibility that the SZ effect can be used as 
a probe of cosmology and only concentrate on the as- 
trophysical uses of the SZ effect. This is a reasonable 
approach to take since there are many unknown astro- 
physics associated with the SZ effect involving the clus- 
tering of gas density and temperature. Such an approach 
allows us not to complicate the parameter measurements 
by adding both astrophysical and cosmological parame- 
ters. Assuming the cosmology will be safely known, we 
now ask the question what additional astrophysical pa- 
rameters one can hope to extract from the SZ effect under 
the present halo model. 

There are many approaches to parameterize the un- 
known astrophysics of the SZ power spectrum. Some 
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possibilities have already been suggested in the literature, 
essentially involving the gas evolution Q]. Since the SZ 
effect involves both gas and temperature as a product, ie. 
the pressure, one may be led to conclude that it is not 
possible to separate effects associated with temperature 
from those associated with gas density. Given the depen- 
dence of temperature on the pressure profile, which we 
calculated from hydrostatic equillibrium, one can expect 
the degeneracy between gas properties and temperature 
properties to be partly broken. Note that the profile 
shape only depends on the temperature and not the ex- 
act density of gas; thus, profile shape is only affected by 
the temperature. This statement is exact under the as- 
sumption that viral equation can be applied. If one were 
to consider a varying temperature distribution, then the 
pressure profile will be sensitive to both gas density and 
temperature. We do not consider the case with a tem- 



perature profile as we are only interested in the most 
simplest case to describe pressure here. 

As discussed earlier, the clustering of pressure power 
spectrum has a turnover corresponding to an equivalent 
scale radius of pressure. Through the gas pressure profile, 
this turnover can be characterized by the parameter b and 
the dark matter scale radius r^. Note that 6 cx l/Tg, so 
a measurement of b is essentially a probe of the electron 
temperature, though, it is unlikely that one can obtain 
all information on temperature and its evolution from 
one parameter measurement. Thus, instead of b, we take 
temperature itself to be one interesting astrophysical pa- 
rameter and consider its evolution such that 

(64) 
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FIG. 7. The non-Gaussian correlation coefficient Cij'^, of 
tlie SZ power spectrum, involving only the configuration of 
the SZ trispectrum that contribute to the SZ power spectrum 
covariance (see, Eq. ^3|). The correlations are such that they 
tend to 1 as — > Ij and is fully described by the contribution 
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to the trispectrum by the single halo term. 



Here, the four parameters represent the temperature- 
mass normalization, Tq, which in the fiducial case has 
the value given by the virial equation (KbTq ~ 5.2 keV), 
the mass dependence slope, Ti, with a fiducial value of 
2/3, a redshift dependent evolutionary parameter, Tcvoi, 
with a fiducial value of 1, and a minimum temperature 
for gas independent of mass and redshift Tmin, with a 
value of zero in the fiducial case. This latter param- 
eter accounts for any possible preheating of gas before 
virializing in halos due to effects associated with some 
unknown astrophysics, such as the reionization process. 
A measurement of Tmin would be interesting given that 
observational data from massive clusters to less-massive 
galaxy groups suggest possible preheating of gas before 
virialization in halos. As discussed in and the 
SZ power spectrum provides a strong probe of preheating 
of gas. We note here that our description of preheating 
as a redshift independent temperature is likely to be too 
extreme since one expects the preheating temperature 
to vary with both redshift and mass such that all three 
parameters, Tq, Ti and Tevoi, are affected. We do not 
consider variations in Tmin as a function of redshift due 
to our unknown knowledge of physics associated with this 
process and that addition of extra parameters can lead 
to further degeneracies. In Fig. |^, we show the variation 
in the redshift evolution of the density weighted temper- 
ature of electrons about the fiducial model. The density 
weight temperature was calculated following Eq. ^ 

In addition to the temperature, the SZ effect also 
depends on the number density of electrons in clus- 



1+z 

FIG. 8. The variation in the density weighted temperature 
of electrons as a function of redshift. The solid line shows the 
redshift evolution of the temperature under the fiducial model 
while variations about this model are shown as labeled. 



ters. So far, we have considered this number through 
the universal baryon fraction in the universe such that 
fg = Mg/Ms = ilg/firn- TMs assumptloH ignores any 
possible effects associated with the evolution of the gas 
fraction in halos, independent of any evolution that may 
be associated with temperature. It will be interesting to 
study to what extent future observations will allow the 
measurement of the fraction of baryons that is respon- 
sible for the SZ effect, and any evolution that may be 
associated with this fraction. Thus, a second set of pa- 
rameters one can hope to extract from SZ observations 
involves gas mass fraction of halos and its evolution. 

To study such gas properties, we parameterize the gas 
mass fraction such that 



fg — fo 



M 







fi 



(i + zy 



(65) 



In a recent paper, Majumdar |Q has suggested the pos- 
sible measurement of any mass and redshift dependence 
of gas mass fraction in galaxy clusters, given that the SZ 
power spectrum was observed to vary significantly with 
changes in these two parameters. Since the SZ power 
spectrum essentially is sensitive to ^ fg^^e^ however, 
such a suggestion for measurement of gas evolution is 
not independent of any variations associated with tem- 
perature, which was ignored in the study of p^ . Our 
general parameterization above involving both tempera- 
ture and gas allows us to quantify how well independent 
statements can be made on possible measurement of gas 
density and temperature evolution, under the assump- 
tion that cosmology is known. Note that gas evolution 
is not present in our fiducial model since we take the 
gas fraction to be independent of mass and redshift with 
/i = and /cvoi = 0, respectively. 
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We now have a total of seven parameters we wish to ex- 
tract from a measurement of the SZ power spectrum. In 
order to perform this calculation we take a Fisher matrix 
based approach. The Fisher matrix is simply a projection 
of the covariance matrix onto the basis of astrophysical 
parameters pi 



TABLE III. 

Inverse Fisher Matrix (x 10^) 



' a/3 



(66) 



where the total covariance includes both the signal and 
noise covariance. Under the approximation of Gaussian 
shot noise, this reduces to replacing Cf ^ Cf^ + Cf°^^° 
in the expressions leading up to the covariance equa- 
tion (|60[). The noise power spectrum includes the noise 
associated with detectors, beam size and variance result- 
ing from the separation of the SZ effect from other tem- 
perature fluctuations in multifrequency data. 

Under the approximation that there are a sufficient 
number of modes in the band powers that the distribu- 
tion of power spectrum estimates is approximately Gaus- 
sian, the Fisher matrix quantifies the best possible errors 
on cosmological parameters that can be achieved by a 
given survey. In particular is the optimal covariance 
matrix of the parameters and {F^^)]^'^ is the optimal er- 
ror on the ith parameter. Implicit in this approximation 
of the Fisher matrix is the neglect of information from 
the parameter dependence of the covariance matrix of 
the band powers themselves. We neglect this informa- 
tion due to computational restrictions on the calculation 
of covariance for all variations in parameters within a 
reasonable amount of time. We do not expect this ex- 
clusion to change our results significantly. Also, here, 
we are mostly interested in an order of magnitude esti- 
mate on how well the SZ power spectrum can constrain 
astrophysics associated with large scale pressure. 

The Fisher matrix approach allows us to address how 
well degeneracies are broken in the parameter space un- 
der the assumption of a fiducial model for the parame- 
ters. For the purpose of this calculation, we take binned 
measurements of the SZ power spectrum following the 
binning scheme in Table 2. We consider a perfect, no 
instrumental noise, SZ experiment with no noise contri- 
bution to the covariance and observations out to I ~ 10^. 
To consider a real world scenario, we also study the as- 
trophysical uses of the SZ power spectrum that can be 
extracted from the Planck mission. Here, we use the SZ 
noise power spectrum calculated for Planck with detec- 
tor noise and uncertainties in the separation of SZ from 
CMB and other foreground in This SZ noise power 
spectrum is shown in Fig. H(b). 



D. Discussion 

In Fig. I, 

we show the variation associated with SZ temperature 
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NOTES. — Inverse Fisher matrix for the SZ effect with seven pa- 
rameters and full non-Gaussian errors. The error on an individual 
parameter is the square root of the diagonal element of the inverse- 
Fisher matrix for the parameter while off-diagonal entries of the 
inverse Fisher matrix shows correlations, and, thus, degeneracies, 
between parameters. We have assumed a perfect, no instrumental 
noise, experiment with a full sky survey (/sky = !)• The seven 
parameters are described in § IV C . 



fluctuations (AT = + 1) / {2tt)CiTcmb) for parame- 
ters involved with gas and temperature evolutions. These 
plots allow us to understand some of the degeneracies as- 
sociated with the description of gas and temperature evo- 
lution. For example, as shown, the gas and temperature 
redshift dependence essentially predicts similar behavior 
for the SZ temperature fluctuation, though there are mi- 
nor differences due to the temperature dependence on the 
pressure profiles of halos. For the most part, variations 
due to temperature evolution is due to the normaliza- 
tion and not due to variations in the profile shape. In 
(b) and (e), we show variations due to the mass slope 
of the gas evolution and temperature evolution, respec- 
tively. Here again, we see similar behavior. When the 
slope of the mass-temperature relation, as a function of 
mass, is greater than 0.7, we see significant differences, es- 
pecially involving an increase in temperature fluctuations 
at small scales. This is due to the relatively increasing 
weighing of small mass halos. 

In (c) and (f), we show variations associated with gas 
evolution normalization /o and temperature-mass nor- 
malization Tq. The variation associated with /o is easily 
understood since the effect is only a change in the overall 
normalization of the power spectrum. The variation with 
temperature-mass normalization shows both effects due 
to normalization and the profile. When the normaliza- 
tion is low, gas clusters to small radii in low mass halos 
leading to an increase in power at small scales. As the 
temperature normalization is increased, gas profile varies 
such that there is a reduction in small scale power and the 
angular multipole of the turn-over scale shifts to low val- 
ues. When the temperature normalization is sufficiently 
high, the overall weighing resulting from the overall tem- 
perature multiplicative factor becomes important. Now, 
the power spectrum behaves as a simple normalization 
change, similar to the variation in power due to gas evo- 
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lution normalization. As shown in Fig. |^(a) to (f), there 
are significant degeneracies involved with astrophysical 
parameters that lead to the SZ effect. 

In Table [V, we tabulate the errors on these seven pa- 
rameters using the inverse Fisher matrix for a possible 
SZ power spectrum measurement. Here, we have con- 
sidered the possibility that parameter extraction will be 
limited to 3, 5 and 7 parameters. The increase in number 
of parameters to be measured from a SZ power spectrum 
increases degeneracies associated with the set of param- 
eters resulting in their accuracies. In the case of the 
3 parameters involving temperature-mass normalization, 
To, a minimum temperature for all halos Tmin, and the 
gas mass fraction /o, in a perfect experiment, all three 
parameters can be extracted such that they will be pro- 
vide essentially very strong constraints. For example, the 
error on /q is such that one can identify the gas fraction 
of clusters responsible for SZ effect from the cosmic mean 
of ng/rim = 0.05/0.3 with an error of 4 x lO'^. With 
Planck, one can constrain the preheating temperature at 
the level of ^ 0.7 keV, and since current predictions for 
possible preheating is also at the level of few tenths keV, 
Planck SZ power spectrum can either confirm or put a 
useful limit on preheating temperature at current expec- 
tations. 

As tabulated, however, the accuracy to which parame- 
ters can be determined from SZ power spectrum reduces 
significantly when the number of parameters to be de- 
termined is increased. For example, the Planck mission 
will only set a limit at ~ 1.7 keV on preheating, if one 
were to study both the mass and redshift dependence of 
electron temperatures. Such an upper limit is unlikely to 
be useful for current studies related to preheating of gas. 
Given that we cannot obtain useful errors with Planck for 
5 parameters, we suggest that Planck may not be useful 
for the purpose of studying the full parameter space sug- 
gested here. This is understandable since Planck only 
allows the measurement of the SZ power spectrum out 
to I 1500, while most of the variations due to pa- 
rameters under discussion here happens at I ~ 5000 or 
higher where the turnover in the SZ power spectrum is 
observed. The Planck mission, however, allows one to 
obtain reasonable errors on parameters which generally 
define the normalization of the power spectrum, such as 
the temperature-mass normalization or the normaliza- 
tion of gas mass fraction. The normalization for gas mass 
fraction from Planck will be useful for the purpose of 
understanding what fraction of cosmic baryons reside in 
massive halos and contribute to the SZ effect and to look 
for any discrepancy of such a value from the total baryon 
content predicted by big bang nucleosynthesis arguments. 
In order to obtain reliable measurements of evolution of 
gas and temperature, a small scale experiment sensitive 
to multipoles out to I ~ 10* will be necessary. 

For a perfect experiment, we show the errors on seven 
parameters also in Table |^ The inverse Fisher matrix 
in this case is tabulated in Table HI. The diagonals of 



parameters, while, more importantly, the off diagonals 
show the covariance between parameters. These covari- 
ances allow one to understand the degeneracies between 
parameters. In Table we show the full extent to which 
parameters degrade the accuracies by tabulating degra- 
dation factors associated with the seven parameters. The 
degradation factor list the increase in parameter error 
from what can be achieved if all other parameters are 
known to what can be achieved when all parameters are 
to be retrieved from data. The degradation factors are 
at the level of one hundred or more for some param- 
eters, suggesting that there are significant degeneracies 
associated with the parameterization of the temperature 
and gas fraction as a function of mass and redshift. Our 
results generally suggest that accurate estimates of gas 
evolutionary properties, in the presence of unknown tem- 
perature properties, are not possible. 

In addition to the parameter degeneracies, the non- 
Gaussianities associated with the SZ effect also increase 
the errors on parameters. For example, for the seven 
parameters under discussion here and again for a per- 
fect and full sky experiment, we list the errors on pa- 
rameters one can obtain if one were to ignore the non- 
Gaussian contributions to the covariance. As tabulated, 
non-Gaussianities increase the error on parameters by up 
to factors of 1.5, suggesting that the ignoring the non- 
Gaussianities will lead to a significant underestimate of 
the errors in parameters. This should be considered un- 
der the context that the SZ effect is significantly non- 
Gaussian at all scales of interest and that ability to dis- 
tinguish parameters happen only at multipoles of a few 
thousand where the non-Gaussianities in fact dominate. 



V. SZ KINETIC (OSTRIKER-VISHNIAC) 
EFFECT 

The bulk flow of electrons that scatter the CMB pho- 
tons, in the reionized epoch, lead to temperature fluctu- 
ations through the well known Doppler effect: 



r^°p(n) 



drg{r)h ■ Vg(r, nr) , 



(67) 



where is the baryon velocity. In Fig. ^(a), we show 
the linear Doppler effect, including contributions result- 
ing due to double scattering effect described in Q (see, 

for details). The power spectrum is such that it peaks 
around the horizon at the reionization projected on the 
sky today. The effect cancels out significantly at scales 
smaller than the horizon at scattering since photons scat- 
ter against the crests and troughs of the perturbation. 

The Ostrikcr-Vishniac effect arises from the second- 
order modulation of the Doppler effect by density fluctu- 
ations 1 22 1, and avoids strong cancelation associated with 
the linear Doppler effect. This nonlinear effect is also 
known as the kinetic SZ effect from large-scale structure 



the inverse Fisher matrix show the variance of individual 
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TABLE IV. 
Parameter Errors 





N 


To 


Ti 




Tmin 


fn 






3 


Perfect 


0.04 






0.002 


0.0004 








Planck 


0.79 






0.75 


0.03 






5 


Perfect 


0.13 


0.02 


0.06 


0.05 


0.002 








Planck 


1.39 


0.41 


1.22 


1.37 


0.05 






7 


Perfect 


0.30 


0.07 


0.17 


0.07 


0.008 


0.06 


0.13 




Degradation 


47 


184 


133 


34 


82 


240 


130 




Gaussian 


0.18 


0.04 


0.10 


0.04 


0.005 


0.04 


0.08 




Increase (%) 


64 


70 


77 


81 


73 


60 


68 



NOTES. — Parameter errors, iF~^)^^'^ , using the halo model and 
the full covariance for the SZ effect. We tabulate these errors for a 
perfect experiment with no instrumental noise and full sky obser- 
vations out to I ~ 10^. We also show the expected errors for Planck 
mission with a useful sky fraction of 65% (/sky = 0.65), and with 
the noise pow er spe ctrum shown in Fig. ^(b). The parameters are 
described in I 
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We break the parameter estimation to consider 
recovery of 3, 5 and 7 parameters. Under "Degradation" we tabu- 



late the degradation factors, {F)ii^^^ /(F~^)'^l^ , due to parameter 
degeneracies. We also list the parameter errors expected if one 
were to assume Gaussian sample variance only for the SZ power 
spectrum and were to ignore the non-Gaussian covariance. The 
increase in error on individual parameters, with the introduction 
of the full covariance matrix, ranges from 40% to nearly 100%. 
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|P3| and is associated with the hne of sight motion of ha- 
los. Due to the density weighting, the kinetic SZ effect 
peaks at small scales: sub arcminutes for ACDM. For 
a fully ionized universe, contributions are broadly dis- 
tributed in redshift so that the power spectra are moder- 
ately dependent on the optical depth r. Here, we assume 
an optical depth to ionization up to 0.1, consistent with 
current upper limits on the reionization redshift from 
CMB [Q and other observational data (see, e.g., | |47| ] 
and references therein). 

The kinetic SZ temperature fluctuations, denoted as 
kSZ, can be written as a product of the line of sight 
velocity, under linear theory, and density fluctuations 

T^^^{h) = J drg{r)n ■ Vg{r, nr)Sg{r, nr) 



drgGG 



(27r)3 

X(5f (k-k')(53(k')e'''-"^ 



n • 



k-k' 
Ik-k 



f\2 



(68) 



Here, we have used the linear theory to obtain the large 
scale velocity field in terms of the linear dark matter 
density field. The multiplication between the velocity 
and density fields in real space has been converted to a 
convolution between the two fields in Fourier space. 

We can now expand out the temperature perturbation 
due to kinetic SZ effect, T'^sz^ jj^^^ spherical harmonics: 



«lf ^ = J dnj dr (gGG) 



d^ki f d^k2 



(27r)3 J (27r)3 
n • ki 



Yrin), (69) 



where we have symmetrizised by using ki and k2 to rep- 
resent k — k' and k' respectively. Using 



n.k = ^|^fcrr'(n)rr'*(k), 

and the Rayleigh expansion 

g.k.fi. ^ 4^^»'ji(fcr)rr*(k)rr(n) 

Im 



(70) 



(71) 



we can further simplify and rewrite the multipole mo- 
ments as 



.-kSZ 



I dr I I S^Y^Hi: 

lirrii l2Tn2 m' 



(27r)3y (27r) 



X {gGG)^-^^^,n, ik,r)S't{^i)SM 

xy™^(ki)rr'(ki)r,r(k2)Cn^"^"''(n) • (^2) 

Here and throughout, we make use of the general integral 
over spherical harmonics written such that 

C™.l-'"'(n)= /rffi>ir^(n)r,r(n)....i^™^(n) (73) 

We can construct the angular power spectrum by con- 
sidering {ai-^rniCLi2m2) ■ Under the assumption that the 
temperature field is statistically isotropic, the correlation 
is independent of m, and we can write the angular power 
spectrum as 



kSZ 



, *,kSZ kSZ \ xD rD ^: 

The correlation can be written using 



(74) 



dngGG / dr2gGG 



,^*,ksz kSz \ _ (47r)^ 

\"iimi "■hm^l — g 

d3ki d3k2 d^k; d^k'^ 
(2^(2^(2^(2^ 

^ (<5f(k;)5,(k^)5f"(ki)5*(k2)) 

X (-^)'l+'l (i)'2+'2j;, L_ j;,,(fc2ri) 

xY/'{k[)Y"''^\k'^)Y;?^k[)Yi^^^^^ 



m2'm2 ■m-2 m2 



/ ^ \ tIIL-, IIL-, IIL-, Ilbf, / ^ \ 



(75) 



We can separate out the contributions such that the to- 
tal is made of correlations following {vgVg){SgSg) and 
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(''^s^ff) (^9*^9) depending on whether we consider cumu- 
lants by combining ki with k'j^ or respectively. Af- 
ter some straightforward but tedious algebra, and noting 
that 



I2 l\ 

^-^ \ m'l m'2 mi 



1'2 



m2 



+ 1 



(76) 



we can write 
22 
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47r 



/ h I2 




drigGG / dr2gGG / fc^dfci / kldh 



X (p]*"(fci)P,,(fc2)j;,(fc2r2)j;,(fc2ri)^ 
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ki 



(77) 



Here, the first term represents the contribution from 
{vgVg){SgSg) whllc thc sccoud term is the {vgSg){vgdg) 
contribution, respectively. In simplifying the integrals 
involving spherical harmonics, we have made use of the 
properties of Clebsh-Gordon coefficients, in particular, 
those involving I — 1. The integral involves two dis- 
tances and two Fourier modes and is summed over the 
Wigner-3j symbol to obtain the power spectrum. The 
above equation represents the angular power spectrum 
of kinetic SZ, or OV effect, under the all-sky corrdinates 
and we have not used the flat-sky assumption or any 
small angular scale limit which is usually considered in 
the literature. 

Since we are primary interested in the contribution 
at small angular scales here, we can ignore the contri- 
bution to the kinetic SZ effect involving the correlation 
between linear density field and baryons and only con- 
sider the contribution that results from baryon-baryon 
and density-density correlations. In fact, under the halo 
description provided here, there is no correlation of the 
baryon field within halos and the velocity field traced 
by individual halos. Thus, contribution to the baryon- 
velocity correlation only comes from the 2-halo term of 
the density field-baryon correlation. This correlation is 
suppressed at small scales and is not a significant con- 
tributor to the kinetic SZ power spectrum (see, [p3| ). 

Similar to the Limber approximation BO], in order 
to simplify the calculation associated with (VgVg) {6 g 6 g), 
we use an equation involving completeness of spherical 
Bessel functions: 



dkk^F{k)3i{kr)ji{kr') « -d'/S'^ir - r')F(k)\.^ 



(78) 



where the assumption is that F{k) is a slowly-varying 
function. Applying this to the integral over fc2 gives 
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The alternative approach, which has been the calcula- 
tional method in many of the previous papers H^,^,^ , 
is to use a specific coordinate frame with the z-axis along 
k. This allows one to simplify the SZ kinetic power spec- 
trum to: 



C, 



kSZ 



1 

8^ 



dr^-^PssikflJk^l- 
di V "-A 



(80) 



with the mode-coupling integral given by 



Ivik) 



dki 



0?/i 



(1 - /i')(l - 2^,yi) Pssjkyi) Ps5iky2: 
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(81) 



We refer the reader to p8[ for details on this deriva- 
tion. In above, = k • ki, j/i = fci//c and y2 = k2/k = 
This flat-sky approximation makes use 



vl 



of the Limber approximation pj| | to further simplify the 
calculation with the replacement of fc = l/dA- The power 
spectra here represent the baryon field power spectrum 
and the velocity field power spectrum; the former as- 
sumed to trace the dark matter density field while the 
latter is generally related to the linear dark matter den- 
sity field through the use of linear theory arguments. 

The correspondence between the flat-sky and all-sky 
formulation can be obtained by noting that in the small 
scale limit contributions to the flat-sky effect comes when 
^2 = |k — ki| ^ k such that j/i ^ 1. In this limit, the 
flat sky Ostriker-Vishniac effect reduces to a simple form 
given by Hu in [^: 



^kSZ 



, (gGG)2 2 

ar ^ggyl^jVm 
"-A 



(82) 



y2 

^r77is 

large scales 



Here, v'i^^ is the rms of the uniform bulk velocity from 



= / dk 



Pssjk) 
27r2 



(83) 



The 1/3 arises from the fact that rms in each component 
is l/3rd of the total velocity. 

In the same small scale limit, to be consistent with the 
flat sky expression, we can reduce the all-sky expression 
such that contributions come from a term that looks like 
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FIG. 10. The tempera- 

ture fluctuation power (AT,^ = l{l + l)/(27r)C;rcMB) for a 
variety of methods to calculate the kinetic SZ effect. Here, 
we show the contribution for a reionization redshift of ~ 8 
and an optical depth to reionization of 0.05. The contribu- 
tions are calculated under the assumption that the baryon 
fleld traces the non-linear dark matter {Pg{k) = P5{k) with 
Ps{k) predicted by the halo model), the linear density field 
(Pg(fc) = P''"(fc)), and the halo model for gas, with total and 
the 2halo contributions shown separately. For the most part, 
the kinetic SZ effect can be described using linear theory, and 
the non-linearities only increase the temperature fluctuation 
power by a factor of a few at / ~ 10^. 
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A comparison of the reduced all-sky (Eq. |7^) and flat- 
sky (Eq. ^) formula in the small-scale limit suggests that 
the correspondence between the two arises when 

^(2Zi + l)(% + l)(o \ \ )\k.{kr)Y (85) 

Numerically, we determined this to be true as long as 
I2 <^ I, however, we have not been able to prove this 
relation analytically. We leave this as a challenge to our 
readers. Note that the limit ^2 ^ ^ is the small-scale limit 
in the flat-sky and is effectively equivalent to j/i ^ 1 
considered in writing Eq. |8^. Thus, the right hand side 
of the above expression denotes the all-sky equivalent of 
the integral that leads to a l/3rd of rms of a randomly 
directed quantity along one particular line of sight. 



A. Discussion 

In Fig. ^, we show our prediction for the SZ kinetic ef- 
fect and a comparison with the SZ thermal contribution. 




FIG. 11. Cumulative signal-to-noise for the detection of SZ 
thermal-SZ thermal-SZ kinetic bispectrum (a) and skewness 
(b) with temperature anisotropy data. The dotted line is 
for the detection of SZ thermal-SZ kinetic correlation using 
CMB data alone, while the dashed line is the same when 
the SZ thermal effect has been separated from other CMB 
contributions and the measurement now involves two points 
from the SZ map and one point from the CMB. Finally, the 
solid line is the maximum signal-to-noise for achievable with 
the separation of the SZ kinetic effect from all contributors to 
CMB anisotropy. 



As shown, the SZ kinetic contribution is roughly an order 
of magnitude smaller than the kinetic SZ contribution. 
There is also a more fundamental difference between the 
two: the SZ thermal effect, due to its dependence on 
highest temperature electrons is more dependent on the 
most massive halos in the universe, while the SZ kinetic 
effect arises more clearly due to large scale correlations of 
the halos that make the large scale structure. The differ- 
ence arises from that fact that kinetic SZ effect is mainly 
due to the baryons and not the temperature weighted 
baryons that trace the pressure responsible for the ther- 
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mal effect. Contributions to the SZ kinetic effect comes 
from baryons tracing all scales and down to small mass 
hales. The difference associated with mass dependence 
between the two effects suggests that a wide- field SZ ther- 
mal effect map and a wide-field SZ kinetic effect map will 
be different from each other in that massive halos, or clus- 
ters, will be clearly visible in a SZ thermal map while the 
large scale structure will be more evident in a SZ kinetic 
effect map. Numerical simulations are in fact consistent 
with this picture [|l6| . 

As shown in Fig^Wb) , the variations in maximum mass 
used in the calculation does not lead to orders of mag- 
nitude changes in the total kinetic SZ contribution; This 
change is less than the changes in the total thermal SZ 
contribution as a function of maximum mass. This again 
is consistent with our basic result that most contributions 
come from the large scale linear velocity modulated by 
baryons in halos. Consequently, while the thermal SZ 
effect is dominated by shot-noise contributions, and is 
heavily affected by the sample variance, the same is not 
true for the kinetic SZ effect. 

In Fig. we show several additional predictions for 
the kinetic SZ effect, following the discussion in [ p3| . 
Here, we have calculated the kinetic SZ power spectrum 
under several assumptions, including the case when gas is 
assumed to trace the non-linear density field and the lin- 
ear density field. We compare predictions based on such 
assumptions to those calculated using the halo model. 
As shown, the halo model calculation shows slightly less 
power than when using the non-linear dark matter den- 
sity field to describe clustering of baryons. This difference 
arises from the fact that baryons do not fully trace the 
dark matter in halos. Due to small differences, one can 
safely use the non-linear dark matter power spectrum to 
describe baryons. Using the linear theory only, however, 
leads to an underestimate of power at a factor of 3 to 4 
at scales corresponding to multipoles of Z ~ 10'' to lO'"^ 
and may not provide an accurate description of the total 
kinetic SZ effect. 

The interesting experimental possibility here is 
whether one can obtain an wide-field map of the SZ ki- 
netic effect. Since it is now well known that the unique 
spectral dependence of the thermal SZ effect can be used 
to separate its contribution [|2j , at smaller angular scales, 
it is likely that after the separation, SZ kinetic effect will 
be the dominant signal, even after accounting for the 
lensed CMB contribution. For such a separation of the 
SZ thermal effect to be carried out and such that a detec- 
tion of the kinetic SZ effect will be possible, observations, 
at multifrequencies, are needed to arcminute scales. Up- 
coming interferometers and similar experiments will al- 
low such studies to be eventually carried out. A wide-field 
kinetic SZ map of the large scale structure will eventu- 
ally allow an understating of the large scale velocity field 
of baryons, as the density fluctuations can be identified 
through cross-correlation of such a map with a similar 
thermal SZ map. We now discuss the existence of corre- 
lations between the SZ thermal and SZ kinetic effect. 



VI. SZ THERMAL-SZ KINETIC CORRELATION 

The SZ thermal and SZ kinetic effects both trace the 
large scale structure baryons. One can study a corre- 
lation between these two effect to probe the manner in 
which baryons are distributed in the large scale struc- 
ture. For example, such a correlation study may allow 
one to answer to what extent diffuse baryons contribute 
to thermal SZ when compared to their contribution to ki- 
netic SZ. Given that the SZ kinetic effect is second order 
in fluctuations, there is no direct two-point correlation 
function between the temperature anisotropics produced 
by SZ thermal and kinetic effects. As discussed in to 
the lowest order, the correlation between kinetic SZ and 
thermal SZ manifests as a nonvanishing bispectrum in 
temperature fluctuations and can be studied by consid- 
ering a three-point correlation function or a bispectrum, 
the Fourier space analog of the three point function, or 
associated statistics such as the third moment or skew- 



A. SZ thermal-SZ thermal- SZ kinetic bispectrum 

We refer the reader to ||^ for full details on the deriva- 
tion of the SZ thermal-SZ kinetic correlation as a bispec- 
trum involving two measurements of SZ thermal effect 
and one measurement of SZ kinetic effect (SZ thermal-SZ 
thermal-SZ kinetic); note that in 1^, we identify this bis- 
pectrum as SZ-SZ-OV. As discussed there, for low optical 
depths to rcionization, the detection of the SZ-SZ-OV 
bispectrum is problematic due to high cosmic variance 
resulting from primary anisotropics. 

In Cooray & Hu @, we assumed that pressure and 
baryons both traces dark matter to calculate the SZ-SZ- 
OV bispectrum. Using our halo model, we can now up- 
date this calculation to include the bispectrum formed 
between SZ thermal-SZ thermal-SZ kinetic effects. Ad- 
ditionally, we can investigate the improvements in the 
signal-to-noise for the bispectrum detections when the 
SZ thermal effect is separated from CMB. The maximum 
signal-to-noise for this bispectrum can only be achieved 
when the SZ kinetic effect is also separated from CMB, 
though, given that the two effects have the same fre- 
quency dependence, it is unlikely that the kinetic SZ ef- 
fect can be separated from thermal CMB anisotropics 
using frequency information alone. Using the frequency 
information, with the thermal SZ contribution separated 
from CMB, however, it is likely that the SZ kinetic effect 
will dominate the small angular scale signal in the tem- 
perature anisotropics. Thus, one can use small angular 
scale thermal CMB temperature fluctuations to construct 
the SZ thermal-SZ kinetic bispectrum. 

Following we can write the SZ thcrmal-SZ thermal- 
SZ kinetic bispectrum as 
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FIG. 12. The baryon-pressure (a) power spectrum (b) 
trispectrum today {z — 0) broken into individual contribu- 
tions under the halo description. The line labeled 'correlation' 
shows the correlation coefficient of gas-pressure correlation 
with respect to gas-gas and pressure-pressure. For reference, 
we also show the pressure power spectrum and the trispec- 
trum. 



with 



^ nli(2^. + i) 



tracing the linear velocity field. Since there is no contri- 
bution to the large scale bulk flows from the non-linear 
regime (ie. the 1-halo term), we model the pressure- 
density cross power as the large scale density-pressure 
correlations in the linear regime described by the 2-halo 
term. 

The signal-to-noise for the detection of the bispectrum 
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Here, "Perm." means a sum over the remaining 5 per- 
mutations of (hjhjh) as usual. 

Here, PgYi is the baryon-pressure power spectrum while 
the Psn is the density-pressure power spectrum, with S 



In |1^, we showed how multifrequency cleaning of SZ 
effect can be a useful tool for higher order correlation 
studies and discussed how the signal-to-noisc for the de- 
tection of SZ-lensing correlation, again through a bispec- 
trum, can be improved by using CMB primary anisotropy 
separated SZ map. We consider the same approach here, 
where we study the possibility for a detection of the 
SZ thermal-SZ kinetic correlation by using a frequency 
cleaned SZ thermal map, which will provide two mea- 
surements, and a CMB map containing only the CMB 
primary, SZ kinetic and and thermal-like secondary con- 
tributions. 

In Fig. ^l](a), we update signal-to-noise results for the 
bispectrum given in ||], where we only studied the pos- 
sible detection in CMB data alone and with no consid- 
eration for separation of effects, especially the SZ ther- 
mal effect. The separation allows a decrease in cosmic 
variance, as the noise is no longer dominated by CMB 
primary anisotropies. This leads to an increase in the 
cumulative signal-to-noise. With SZ thermal effect sepa- 
rated, we see that the signal-to-noise increases by roughly 
two orders of magnitude. In |p^ , we showed how one can 
obtain an order of magnitude in signal-to-noise when a 
CMB separated SZ thermal map is used for a detection of 
SZ-lensing correlation. Here, we obtain roughly two or- 
ders of magnitude improvement, since the SZ thermal-SZ 
kinetic correlation is present with two SZ thermal mea- 
surements, instead of one in the case SZ thermal-lensing 
correlation through the SZ thermal-CMB-CMB bispec- 
trum. 

If one can separate the SZ kinetic such that a perfect 
SZ kinetic map, as well as a perfect SZ thermal map, 
is available, then one can improve the signal-to-noise for 
detection significantly such that a detection is possible. 
Since the SZ kinetic effect is expected to dominate tem- 
perature anisotropies at small angular scales, when SZ 
thermal is removed, an opportunity to detect the SZ 
thermal-SZ kinetic correlation will likely come from small 
angular scale multifrequency experiments. One can also 
improve the possibility of detecting this correlation by 
noting that the configuration for the bispectrum is such 
that it peaks for highly flattened triangles (see, 13]). It is 
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likely that progress in experimental studies will continue 
to a level where such studies will eventually be possible. 




FIG. 13. (a) The SZ thermal-SZ kinetic power spectrum 
of squared temperatures. Here, we show the contribution to 
the power spectrum when only Gaussian terms, ie. power 
spectra of pressure-density correlation, are considered (dot- 
ted line) and when when non-Gaussianities are introduced 
through the pressure- density trispectrum (solid line). In (b), 
we show the cumulative signal-to-noise for the detection of 
the SZ thermal-SZ kinetic squared temperature power spec- 
trum using information in multipoles from 2000 to 10000 and 
assuming no instrumental or any other noise contributions to 
the covariance. The signal-to-noise is calculated assuming the 
power spectrum is measured in CMB data (dot-dashed line), 
with a perfect frequency cleaned SZ thermal map (solid line) 
and with a perfect SZ thermal and SZ kinetic effect maps 
(dotted line). In a dashed line, we show the maximum sig- 
nal-to-noise achievable for the power spectrum, with only a 
Gaussian contribution to the covariance. 

Since the bispectrum may be hard to measure directly 
from observational data, we also consider a real space 
statistic that probes the non-Gaussian information at the 
three point level and consider the third moment (see, p^ ] 
for details). In Fig. ^(b), we show the signal to noise for 
the detection of the third moment. Here, we use a top-hat 



window in multipole space out to Zmax so that direct com- 
parison is possible with the signal-to-noise calculation in- 
volving the bispectrum. As shown, we find that there is 
less signal-to-noise in the skewness when compared to the 
full bispectrum. This results from the fact that bispec- 
trum contains all information at the three point level, 
while a measurement of the third moment leads to a loss 
of information. This can also be understood by noting 
that the signal-to-noise for the bispectrum and skewness 
is such that in the case of the bispectrum signal-to-noise 
is calculated for each mode and summed up while for the 
skewness signal-to-noise is calculated after summing the 
signal and noise separately over all modes. 



B. The SZ Thermar-SZ kinetic^ Power Spectrum 

In addition to the SZ thermal-SZ kinetic-SZ kinetic 
bispectrum, we can introduce higher order correlations 
involving the SZ thermal and SZ kinetic effect that probe 
the correlation between the two. One such a possibility is 
the trispectrum formed by the SZ thermal and SZ kinetic 
effect. Given that we do not have a reliable method to 
measure the bispectrum even, the measurement of such 
a higher order correlations in experimental data is likely 
to be challenging. 

Here, we focus on a statistic that captures the corre- 
lation information coming from higher order, essentially 
from a trispectrum, but is easily measurable in exper- 
imental data since it only involves only a power spec- 
trum. Such a possibility involves the power spectrum of 
squared temperatures instead of the usual temperature 
itself. Our motivation for such a statistic came when 
we inspected the published maps of the large scale SZ 
thermal and SZ kinetic effects in simulations by and 
realized that there is a significant correlation between 
the two effects. Since the temperature fluctuations pro- 
duced by the SZ kinetic effect oscillate between positive 
and negative values depending on the direction of the 
velocity field along the line of sight, as stated earlier, a 
direct two point correlation involving the temperature re- 
sults in no contribution. A non-zero correlation between 
the SZ thermal and SZ kinetic effects still manifests if 
the absolute value of the temperature fluctuation due to 
kinetic SZ effect is considered. Since absolute value of 
temperature is equivalent to squaring the temperature, 
we consider the cross-correlation of SZ thermal and SZ 
kinetic effects involving the power spectrum of squared 
temperatures here. 

In order to calculate the SZ thermaP-SZ kinetic^ power 
spectrum, we can take the all-sky approach presented to 
describe the SZ kinetic effect. The end result, however, is 
numerically cumbersome since it involves 5 dimensional 
sum over a Wigner-6j symbol (see, p9| for a derivation). 
Here, we take a flat-sky approach and derive the thermal 
and kinetic SZ squared power spectrum in the limit that 
the velocity field of the kinetic effect is independent of 
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the baryon fluctuations. 
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FIG. 14. The cumulative signal-to-noise for the detection 
of the thermal SZ-kinetic SZ squared temperature power spec- 
trum. In (a), we consider a large angular scale experiment, 
consistent with Planck, and use multipole information from I 
of 100 to 1500. The cumulative signal-to-noise, even with a 
perfectly cleaned SZ map in this multipole range, is signifi- 
cantly less than 1, suggesting that a detection is not possible. 
In (b), we show the cumulative signal-to-noise for a small 
angular scale experiment with information in the multipole 
range of 2000 to 5000. There is adequate signal-to-noise for 
the detection of the squared temperature power spectrum at 
such small scales even if the SZ effect is not completely sep- 
arated from the thermal CMB contribution. For comparison, 
in dotted lines, we show the signal-to-noise achievable if the 
kinetic SZ effect is separated from CMB, in addition to a SZ 
separated map. 

Following our previous definitions, we define the flat 
sky temperature squared power spectrum as 



(e2kSZ(i)02SZ(l,)^ = i2Tr)HD{y + S^-^^^ , (90) 

where the Fourier transform of the squared temperature 
can be written as a convolution of the temperature trans- 
forms 
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Here, it should be understood that 0^(1) refers to the 
Fourier transform of the square of the temperature rather 
than square of the Fourier transform of temperature, 
[8(1)]^. To compute the square of the SZ thermal and 
SZ kinetic temperature power spectrum, we take 

(e2kSZ(i)02SZ(j/)^ = (2^)25^(1 + l')Ci 



(92) 
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Note that the Fourier transform of the temperature fluc- 
tuations in the flat sky is 



6(1) = j (fee- T{e) . 



(93) 



In the small scale limit, where the density field is sep- 
arated from the velocity field, we can separate the above 
cumulant to two parts with one involving just power spec- 
tra of pressure-baryon cross-correlation and another with 
the trispectrum of pressure and baryon fields. We first 
write down the Gaussian-like piece as 
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where we have taken {{9 ■ v){9' ■ v')) ^ l/iv'^^^ with 1/3 
coming from the fact that only a third of the velocity 
component contribute to the line of sight rms. We can 
now introduce power spectra in above correlators such 
that 

(ekSZ(i^)QkSZ(j _ l^)eSZ(i^)0SZ(j, „ i^)^G ^ 
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The integrals over the hne-of-sight wavevectors behave 
such that only perpendicular Fourier modes contribute 
to the projected field, such that l'^ /d\ 3> k"^ ■ This is 
the so-called Limber approximation |Q . Doing the inte- 
gral over the wavevector, then, results in a delta function 
in (ri — such that only contributions come from the 
same rcdshift. Putting the correlator back in the power 
spectrum equation (Eq. |9^), we now get 
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where we have introduced a factor of 2 account for the 
additional permutation involved in the baryon density- 
pressure correlation. 

Similarly, the non-Gaussian piece follows as 
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Since we only use the single halo term to calculate TgUgn, 
the arguments are simply scalars and does not dependent 
on the orientation of the quadrilateral. In general, how- 
ever, the trispectrum depends on the length of the four 
sides plus the orientation of at least one of the diagonals. 



1. Signal-to-Noise 

In order to calculate the possibility for a detection of 
the thermal SZ^-kinetic SZ^ power spectrum, we need 
the covariance of the estimator involved with the mea- 
surement of the squared power spectrum: 
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Here, As = J d^l is the area in the 2D shell in Fourier 
space over which the integral is done and Af is the total 
area of the survey in Fourier space and can be written as 
Af — (27r)2/r2 with a total survey area on the sky of SI. 
Following 1^, we can write down the covariance of our 
estimator as 
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where Cf^^'^ ^^'^ is the squared power spectrum of ki- 
netic SZ and thermal SZ while cf^sz-ksz (^2SZ-sz 
are the squared power spectra of kinetic SZ and ther- 
mal SZ respectively. Here, we assume that squared fields 



are Gaussian. To calculate C, 
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make several assumptions: We assume that the temper- 
ature squared power spectrum will be measured using 
two maps involving frequency separated SZ contribution 
(which will have + 0"°""=) and a map with kinetic 
SZ contribution with the CMB primary component, such 
that the will be composed of eP"'""^ + 0^^^ + . 
Following such a separation, we can write the Cf 
as 
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and like wise for (^zksz-ksz jjgj.g^ contributions come 
from a Gaussian part involving SZ power spectra and 
a non-Gaussian part through the SZ trispectrum. Since 
the primary component fiuctuations dominate the kinetic 
SZ temperature, and that there is no measurable trispec- 
trum for this component under current adiabatic GDM 
predictions we ignore any non-Gaussian contribution to 
Cj and write it as the one with the Gaussian 

part only. This assumption is also safe at small angu- 
lar scales out to I ~ 10^, under the halo model, when 
Qksz > Qprimary ^^^^^ ^j^g kiuctic SZ cffcct Can be de- 
scribed for the most part using large scale correlations 
instead of the single halo term. 

the ratio of As/Af is the total num- 



In the Eq. 99 



ber of modes that measures the squared power spec- 
trum independently and can be approximated such that 
As/Af — /sky (2^ + 1). To calculate the signal-to-noise in- 
volved in the detection of the squared temperature power 
spectrum, we consider an optimized estimator with a 
weighing factor Wi such that 
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and write the signal-to-noise as 
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The weight Wi that maximizes the signal-to-noise is 

the required signal-to-noise as 
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C. Discussion 

In Fig. p^(a), we show the power spectrum of squared 
temperatures for the SZ thermal and SZ kinetic effects 
using the halo term. Here, we have separated the Gaus- 
sian and non-Gaussian contribution to the squared power 
spectrum. As shown, the non-Gaussian contribution to 
the power spectrum is significantly higher than the Gaus- 
sian contributions. 

The Gaussian contribution to the squared power spec- 
trum traces the pressure-baryon density field power spec- 
trum, which is shown in Fig. 12(a) using the halo model. 
In the same figure, for comparison, we also show the 
pressure-pressure power spectrum and the correlation co- 
efficient for pressure-baryon with respect to the pressure- 
pressure and baryon-baryon power spectra. The cor- 
relation behaves such that pressure and baryons trace 
each other at very large scales while the correlation is 
decreased at small scales due to the turnover in the pres- 
sure power spectrum. This is equivalent to the statement 
that there is no low mass halo contribution to the pres- 
sure power spectrum; these halos continue to contribute 
to the baryon density field power spectrum. 

The non-Gaussian contribution to the thermal SZ- 
kinetic SZ squared temperature power spectrum traces 
the trispectrum formed by pressure and density field. 
We show this in Fig. |l2(b) following the halo model. 
For comparison, we also show the trispectrum formed by 
pressure alone in the same figure. The pressure-baryon 
trispectrum is such that at large scales, corresponding to 
linear scales, significant contributions come from the cor- 
relations between halos instead of the single halo term. If 
there are significant contributions coming to the squared 
temperature power spectrum from such linear scales, the 
Gaussian part of the power spectrum should dominate. 
Since all contributions to the squared temperature power 
spectrum comes from small angular scales corresponding 
to non-linear scales in the pressure-baryon trispectrum, 
we only use the single halo contribution in calculating 
the non-Gaussian part of the squared temperature power 
spectrum. In both Gaussian and non-Gaussian parts of 
the power spectrum, the velocity field of the halos are 
taken to be the large scale bulk flows through the linear 
theory. 



In order to assess the maximum possibility for a mea- 
surement of the temperature squared power spectrum in- 
volving kinetic SZ and thermal SZ, here, we ignore the 
detector and beam noise contributions to the covariance. 
Also, we assume full -sky experiments with /sky = 1. 



As written in Eq. 100 



the contribution to the covari- 
ance comes as a convolution in Fourier space. Thus, even 
at small angular scales corresponding to high multipoles, 
noise contributions can come from large angular scales 
or low wavelength modes. Such modes do not have any 
signal and by only contributing to the variance, they can 
reduce the effective signal-to-noise in the measurement. 
Since the squared power spectrum effectively peaks at 
multipoles of ^ 10^, we can essentially ignore any con- 
tribution to the signal, as well as noise, from multipoles 
less than few thousand. These are the same multipoles in 
which the CMB primary anisotropics dominate, thereby, 
increasing the effective noise for the measurement. In 
order to remove the low multipoles, we introduce a filter- 
ing scheme to the spherical or Fourier transform of the 
temperature anisotropy measurements and suppress the 
low multipole data such that the filter essentially acts as 
a high pass filter above some I > /min- 

In Fig. |l3|(b), we show the cumulative signal to 
noise for the measurement of the thermal SZ-kinetic SZ 
squared temperature power spectrum. Here, we have as- 
sumed an experiment, with no instrumental noise, such 
that information is only used in multipoles of 2000 to 
10000. We find significant signal-to-noise for the detec- 
tion of the squared temperature power spectrum, espe- 
cially when using a frequency cleaned SZ map with a 
CMB map, which has no SZ contribution. One can in 
fact use the CMB map itself especially if multifrequency 
information is not available for SZ separation. Since an 
experiment with multipolar information out to I ^ 10'* 
will not readily be available, we consider two separated 
realistic cases, involving a large angular scale experiment, 
similar to Planck, and a small angular scale experiment 
similar to the ones proposed for the study of SZ effect. 
We summarize our results in Fig. ^ As shown in (a), an 
experiment only sensitive to multipolar information rang- 
ing from 100 to 1500 does not have any signal-to-noise 
for a detection of the squared power spectrum. With a 
perfect SZ separated map in this multipolar range, the 
cumulative signal-to-noise for the squared power spec- 
trum is in the order of 0.01. Thus, it is unlikely that 
Planck data will be useful for this study. Since we have 
not included any instrumental noise in calculating signal- 
to-noise, the realistic signal-to-noise for Planck would be 
even lower. 

Going to smaller angular scales, we find that the signal 
increases significantly such that an experiment only sensi- 
tive to the range of I ~ 2000 to 5000 has adequate signal- 
to-noise for a detection of the squared power spectrum. 
A multifrequency experiment in the arcminute scales can 
use its frequency cleaned SZ map to cross-correlate with 
the CMB map and obtain the squared power spectrum 
with a cumulative signal-to- noise of order few tens. A 
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comparison to Fig. |T^(b), suggests that going to lower 
scales beyond 5000 increases the signal-to-noise, and this 
is due to the fact that SZ thermal-SZ kinetic squared 
power spectrum peaks at multipoles of ^ 7000 to 8000, 
suggesting that for an optimal detection of the squared 
power spectrum, one should also include observations out 
to such high multipoles. 

VII. SUMMARY & CONCLUSIONS 

We have discussed non-Gaussian effects associated 
with local large-scale structure contributions to the 
Cosmic Microwave Background temperature fluctua- 
tions during the reionized epoch involving the thermal 
Sunyaev-Zel'dovich effect. At low redshift, distribution 
functions associated with the large scale structure prop- 
erties are non-Gaussian, primarily due to the non-linear 
gravitational evolution. The non-Gaussianities associ- 
ated with the SZ effect is due to the non-Gaussian distri- 
bution of the large scale pressure fluctuations. Since we 
do not have reliable analytical methods to calculate the 
higher order clustering of pressure, or as a matter of fact, 
any other property associated with large scale structure 
including the dark matter distribution, we have utilized 
the so-called halo approach to clustering. The basic de- 
scription of this technique is that correlation functions, 
and their Fourier analogies, can be described through 
clustering of the property within and between halos. The 
halos themselves are assumed to be clustered with respect 
to the linear density field and with a mass dependent bias. 
We use linear theory and its perturbations out to second 
order to describe the clustering of halos. We use the halo 
method to describe clustering of pressure assuming that 
gas is in hydrostatic equilibrium with the dark matter 
density profile in halos. We use the virial equation to 
describe the electron temperatures. 

It is well know that the frequency dependence of the SZ 
effect allows a separation of its contribution from other 
temperature fluctuations and foreground anisotropics in 
multifrequency CMB experiments. Using the pressure 
trispectrum, the Fourier analogue of the four-point cor- 
relation function, calculated under the halo model, we 
discuss the full covariance of the SZ thermal power 
spectrum. The covariance, with the inclusion of non- 
Gaussianities, allows us to properly establish the errors 
associated with SZ power spectrum measurements. Thus, 
we can use this full covariance to study the astrophysical 
uses of the SZ effect. Through a Fisher matrix formal- 
ism, we discuss how well a measurement of the SZ power 
spectrum can be used to conclude details on gas and tem- 
perature evolution and show that there are significant 
degeneracies between gas and temperature evolution pa- 
rameters. The SZ power spectrum from Planck mission 
can be used to understand the amount of baryons present 
in halos that contribute to the SZ effect and to establish 
any nongravitational heating of electrons, such as due to 



preheating with an error of 0.7 keV. 

With the SZ thermal effect separated in temperature 
fluctuations using its frequency information, the kinetic 
SZ effect, also known as the Ostriker-Vishniac, is ex- 
pected to dominate the thermal contribution at small 
angular scales corresponding to multipoles of few thou- 
sand. This effect arises from the baryon modulation of 
the first order Doppler effect. The presence of the SZ ki- 
netic effect can be determined through a cross-correlation 
between the SZ thermal and a CMB map at small scales. 
Since the SZ kinetic effect is second order, however, con- 
tributions to such a cross-correlation arise to the lower 
order in the form of a three-point correlation function, 
or a bispectrum in Fourier space. We suggest an addi- 
tional statistic that can be used to study the correlation 
between pressure traced by the SZ thermal effect and 
the baryons traced by the SZ kinetic effect involving the 
temperature anisotropy power spectrum of squared tem- 
peratures instead of the usual temperature itself. This 
power spectrum probes the trispectrum formed by the 
pressure-baryon cross-correlation. Through a signal-to- 
noise calculation, we show that future small angular scale 
multi-frequency CMB experiments, sensitive to multi- 
poles of a few thousand, will be able to measure the cross- 
correlation of SZ thermal and SZ kinetic effect through 
a temperature squared power spectrum. 
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